


DEC 14°45 NT) Cope 


THE QUARTERLY JOURNAL OF 


MATHEMATICS 


OXFORD SERIES 





Volume13 Nos.50-1 June-Sept. 194.2 





CONTENTS 
T. M. MacRobert: Some Integrals involving 
E-functions and Confluent Hypergeometric 
Functions . ° , . ° . 65 
F. H. Jackson: On Basic Double Hypergeometric 
Functions ‘ ; ; 2 ° . 69 
R. Rado: A Theorem on Independence Relations . 83 
J.L.Burchnall : Differential Equations associated with 
Hypergeometric Functions , : . go 


A. Erdélyi: Integral Equations for Heun Functions . 107 


T. Vijayaraghavan: On Functions represented by 
Certain Series . . ‘ ‘ : 113 


H. L. Hamburger: Contributions to the Theory of 
Closed Hermitian Transformations of Defi- 
ciency-Index (m, m) . ; ; é ° 117 





OXFORD 
AT THE CLARENDON PRESS 
1942 


Price 15s. net 


PRINTED IN GREAT BRITAIN BY JOHN JOHNSON AT THE OXFORD UNIVERSITY PRESS 











THE QUARTERLY JOURNAL OF 
MATHEMATICS 


OXFORD SERIES 


Edited by T. W. CHAUNDY, U. S. HASLAM-JONES, 
J. H. C. THOMPSON 
With the co-operation of A. L. DIXON, W. L. FERRAR, G. H. HARDY, 
E. A. MILNE, E. C. TITCHMARSH 


HE QUARTERLY JOURNAL OF MATHEMATICS 
(OXFORD SERIES) is published in March, June, 
September, and December, at u price of 7s. 6d. net for a single 
number with an annual subscription (for four numbers) of 
278s. 6d. post free. 
Papers, of a length normally not exceeding 20 printed pages 
of the Journal, are invited on subjects of Pure and Applied 


Mathematics, and should be addressed ‘ The Editors, Quarterly 
Journal of Mathematics, Clarendon Press, Oxford’. Con- 
tributions can be accepted in French and German, if in 
typescript (formulae excepted). While every care is taken of 
manuscripts submitted for publication, the Publisher and the 
Editors cannot hold themselves responsible for any loss or 
damage. Authors are advised during the emergency to retain 
a copy of anything they may send for publication. Authors of 
papers printed in the Quarterly Journal will be entitled to 
50 free offprints. Correspondence on the subject-matter of the 
Quarterly Journal should be addressed, as above, to ‘The 
Editors’, at the Clarendon Press. All other correspondence 
should be addressed to the Publisher (Humphrey Milford, 
Oxford University Press, at the temporary address). 


HUMPHREY MILFORD 
OXFORD UNIVERSITY PRESS 
AMEN HOUSE, LONDON, E.C.4 


Temporary address : 
SOUTHFIELD HOUSE, HILL TOP ROAD, OXFORD 





SOME INTEGRALS INVOLVING £-FUNCTIONS 
AND CONFLUENT HYPERGEOMETRIC 
FUNCTIONS 


By T. M. MACROBERT (Glasgow) 


[Received 8 September 1941] 


1. Introductory 

In this paper some integrals involving products of simple cases of 
E-functions and related integrals of products of confluent hyper- 
geometric functions are evaluated. In §2 a list of known formulae 
required in the proofs is given. An integral involving one £-function 
is discussed in §3; and the integrals involving products of the func- 
tions are then derived in § 4. 


2. List of formulae 

The £-function was defined in the Proceedings of the Royal Society 
of Edinburgh, 58 (1937), 3. Only simple cases of the function are 
employed in this paper, the following four formulae alone being 


required: 


@ 


E(a,B::2) = ria) { enea(1 42) dA (R(f) > 0), (1) 
: 
0 


B(a,B 352) = STB—a)P (abe (ai a—B+ 152), (2) 


E(4—k—m, $—k+m:: x) 
= P(j—k—m)P(}—k-+-m)e*el"W, (2), (3) 
cos nw E(4-+-n, 4—n: : 2x) = 4/(27rx)e*K,, (x). (4) 
Note. In (2) the notation >} indicates that to the expression fol- 


a,B 
lowing the symbol is to be added a similar expression with « and 8 


interchanged. 
In addition, use will be made of the well-known formula for the 
hypergeometric function 
, P(y—a—B)l(y) ( a, 8B; ) 
F(a, B3y3;x) = ; ——_F 1—a2)+ 
(oBi7i2) = ya y—p) lat P—yt1 
P(a+fp—y)I(y) _«-pplY—% ¥—B; 
-(l—2x)”Y ss BF 1—x ? (5) 
+ Ta) y—o—B+1 


3695.13 F 
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and of Dixon’s theorem 


1 bs b, = 1 
“|l+a—b, 1+a—c 
C(1+4a)C(1+a—6)(1+a—c)(1+4a—b—c) (6) 
~ Pd+a)P(1+4a—6) (1+ $a—c)P(1+a—b—e)’ 


where R(2+-a—2b—2c) > 0. 
3. An integral involving one #-function 
If l= [ e—t”v-1 (a, B: : ta) dt, 
0 
where x > 0, R(a+y) > 0, R(8+y) > 0, on substituting for the 


E-function from formula (1) and then putting A = pt and changing 
the order of integration it is found that 


I = I'(a) [we a(a+#) dy | et +wyb+y-1 dt 
0 


@ 


_ Ter pet i" Ip. 
( ) Bn | ay (1+ p)a ap 


Hence, if R(x) > $, R(a+y) > 0, R(8+y) > 0, 
pe P(a)P(B)(a-+y)T(B-+y) ( a, fp; 2-1 
etv-1B (a, B:: tx) dt = ——— F 
atl P(a+B+7) atB+y 2 
(7) 
This is equivalent to the formula 
[ ea-2y'—-1f,,(t) dt 
0 
T(l-+-m+4)l(l—m+}4) ,/l+-m+}4, l—m-+4; 
_- 2 2) F : 2x) (8 
Pd—k+1) ( I—k+1 ) (8) 


given by Goldstein.* To see this, replace 2, t, «, 8, y in (7) by 1/(1—2), 
(l—ax)t, 3—k—m, }—k-+m, k+l respectively, and apply formula (3). 


° 


A proof has also been given by Erdélyi.+ 


* Goldstein, Proc. London Math. Soc. (2) 34 (1932), 114. 
t+ A. Erdélyi, Quart. J. of Math. (Oxford), 10 (1939), 189. 
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4. Integrals involving products of Z-functions 


On applying formula (2) to the first of the H-functions in the 
integral below and then integrating term by term, using formula (7), 
it is found that 


| eV B(a,B:: yt)E : zt) dt 


- 2 al (A) a a (atr)C (at+y+A+r)r (a+y+pt+r) 


sina) r!P(a—B+1+r)P(a+y+aAtp+r) * 





A, B; =) 
x a+r ? el 9) 
' Layee z ( 
where ly|<1, R(z)>4, R(aty+A)>0, R(8+y+A) > 0, 
R(at+y+H) > 0, R(B+y+p) > 0. 
Now, when y > 1, this becomes 


= 9) 


[ ety-1H(a,B::t)E(A, pw: : 2t) dt 


> rT (a+y+h+r) 
"_ 2, mat ee -"' ee ag 


* | earn 
0 


_ aes ae 
= PatB+y+n) y+p)T(B+y+E) x 


ay ee Po) 
[arene 
0 





; > P(B—a)P(a+B+y+p) 1 =F” atytp; 1 ) 


TB +y+w 0B) (1+)? \ a—B+1 1+e 


But, by (5), the expression in the last line is equal to 


1 ms atytpy; v ) 
(1+-v)* 


at B+ytp l+e 
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Hence, on making this substitution, expanding the hypergeometric 
series, and integrating term by term, we have 
e-tY-1E (a, B: : t)E(A, w:: 2t) dt 
0 


—. PMP BLAM a+ytAratytwl bry +H) 


D(o+B+yt+p)P(a+y+A+p) 
> (x),(H)-(a+y+ 1), F( A, w+T; =) 


SL, Vat B+ytmlatytAtm), \atytatutr 2 


(10) 
where R(z) > $, R(a+y+A) > 0, R(8+y+A) > 0, R(at+y+p) > 0, 
R(B+y+ pe) > 0. 

By substituting from (3) and (4) in (9) and (10), integrals involving 
products of confluent hypergeometric functions and of modified 
Bessel functions of the second kind can be obtained. 

For example, in (10) put z = 1 and replace a, B, A, wp by $+k+m, 
$-+-k—m, }—k—m, 4—k-+mrespectively. Then, if R(y+2m-+1) > 0, 


(72H (t)W jem (t) at 


0 


P(y+1)P+2m- yy 1) 
C(3+k+m-+y)C(3—k+m-+y) — 
g: sR Poggi ap h+m, Sheil teal 1]. 
$—k+m-+y, §$+k+m+y 
But, by (6), the hypergeometric function is equal to 


C(3+m-+4y)l(3—k+m-+y)P(3+k+m+y)C(4—m-+}y) 
3 by) T(3 y)T(3 vty) '(4—m+4y 


D(2+2m+y)P(1—k+ ty) PL+k+ dy) Py +1) 


Therefore, if R(y+-2m+1) > 0, 


20 ee Speen V1 mt1 
C(y+1)P(4+m-+}y)C(3 m+ sy) (11) 


| 7H, n(t)Wjzm(t) at = 


0 


20 (1+-k+4y)P(1—k+4y) 





ON BASIC DOUBLE HYPERGEOMETRIC 
FUNCTIONS 


By F. H. JACKSON (Eastbourne) 


[Received 14 January 1942: in revised form 29 April 1942] 


1. Introduction 

[N recent papers by Burchnall and Chaundy* many most interesting 
and novel expansions of Appell’s double hypergeometric functions 
are obtained and discussed. The expansions in question were ob- 
tained by using certain operative-functions, formed from hyper- 
geometric functions in which the parameters were replaced by the 
partial-differential operators x@/éx, yé/dy. Mr. Chaundy suggested 
to the present writer that possibly a similar procedure using basic 
operators might give interesting expansions of new basic double 
hypergeometric functions. This has proved to be the case. As such 
functions have not been studied previously I give in this paper the 
q’, g? partial differential equations satisfied by the basic functions. 
These basic equations reduce in the limiting cases to the ordinary 
linear partial differential equations satisfied by Appell’s functions. 
[ do not consider it necessary to give detailed analysis in the case 
of every one of the many expansions which appear: in one or two 
typical cases only will there be detailed analysis, since the method 
and principles used are the same throughout the paper. 

Basic functions are of increasing interest both in number theory 
and in combinatory analysis. This may be some justification for 
writing the present paper, which throughout is due to the suggestions 
and influence of Burchnall’s and Chaundy’s work. I wish here to 
acknowledge the kind help of Dr. W. N. Bailey in respect of Lemmas 
(24) and (25). 


The Normal and Abnormal Functions 

In the case of basic functions two different forms appear, which 
[ term (i) Normal, (ii) Abnormal, respectively. The abnormality 
arises from the presence of a quadratic solitary factor g” in all the 
terms of the series. The two kinds of series, though very similar in 


* J. L. Burchnall and T. W. Chaundy, Quart. J. of Math. (Oxford), 11 
(1940), 249-70. 
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appearance, differ widely owing to the presence of the quadratic 
solitary factor, which factor arises naturally, and, as will be seen, 
cannot be evaded. 


2. The four normal functions 


The functions 


, . = m+n b m(O" n U/ if 
Oa; b,b’;c; 2, y] = > eat ae ) amy”, 
) m+n 


m=0 n= 


(C)m(C’)n 


2 , l , 
O° a D0 2m c’: a, y| = S a (in) Jan -n( b) m( 0 dn amy”, 


m= = n= 


O° a, a’; b,b';¢;2,y] = > Chale ae In gmy, 


m=0 n=0 (m)! (n)! (C)m+n 
; , (4) m+n()m+n 
(m)! (n)! (c),,(€ a 


m=0 n=0 


in which (a),. a4 (1—g*)(1—g**?)...(l1—g**"-), 


. a. 
Oa, b;c,c’32,y] = amy”, 


are basic analogues of functions discussed by Appell and Kampé de 
Fériet.* 


3. The abnormal functions 
The four abnormal functions differ from the four normal functions 


in the fact that ynginn—)) 


replaces y” in the general term and so in all terms of the series. The 
notation for the abnormal ®® will suffice to show the distinction. 


. Pene ° aes . 7 (2) min(®)m(O’)n oan 
OMa;6,0se3295a] = DD nil mil Clmon 


m=0 n=0 
The abnormality is indicated in the expression on the left by the 
presence of g among the elements. 

In the limiting case when |qg| =1 the series 0”, O® converge 
absolutely when |x|, |y| <1, ®® when |2|+|y|<1, ®® when 
\va|+-|vy| <1. It follows that for the same conditions there is even 
sharper convergence when |g| < 1, which condition will be assumed 


yrgtnin—), 


* P. Appell and J. Kampé de Fériet, Fonctions  hypergéomdériques et hyper- 
sphériques (Paris, 1926). 
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throughout the paper. It is obvious that in the case of the abnormal 
series the presence of the solitary quadratic factor in every term will 
ensure a very rapid and wide range of convergence. 


4. Notation 
As in previous papers by the writer, 


(9), = (1—g?)(1—g**?)...(1—g? +") (6) 
($)n = (1—g*)(1—g#+)...(1—gh ") 


-where 0, 6 = xé/éx, yd/ey, q? = exp(@logq), g¢ = exp(¢logq), 
{G(h)}* = lim(1—g")(1—q"*?)...(l—q"**)_ (|q| <1). (7) 


This function {G(h)}— is not the basic gamma function,* but it serves 
to replace the ordinary gamma function in the basic formulae as 
used in the following analysis: 


7, (h) = SNGM+-0+8) G(h+ 0)Gh+4) 

” G(h+8)G(h+¢) G(h)G(h+0+ 4) 
The notation of A and its inverse V due to Burchnall and Chaundy 
is of great convenience. It is by means of these two operators, 
and their equivalents in series, that the expansions of the basic 
double hypergeometric functions, which follow, are obtained. The 


question arises at once, Can such equivalences as 


A,(h) = (8) 


— (—9)(—$)r Bigs 
<= ~ (r)! (h), or (9) 


be justified? The answer is, Certainly they can. It is sufficient to 
note that the operations on a#”y" by V,(h), A,(h), V,(h)A,(k), respec- 
tively and the operations by the respective series equivalent to the 
A operators, always produce terminating series and finite products, 
for example identities such as . 


V,(h)amy” = ()min amy = p (= MAT Mr myn gedhimsn), (10) 


(h) m(h)n (r)! (A), 


which is the basic form of Vandermonde’s identity. Moreover, the 
infinite expansions obtained are absolutely convergent within the 
ranges stated for the variables x, y, and the base q. 


* F. H. Jackson, Proc. Royal Soc. A, 76 (1905), 127-44; J. E. Littlewood, 
Proc. London Math. Soc. (2) 5 (1907), 395. 
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5. The fundamental transformations for the normal functions 
Since V,(h) transforms (h),,(h),, xy” into (h)min vy", We can write 
Oa; b,b’;c,c’; x,y] = K tt b;c;x)@(a, b’;c’; y) 
D®[a, a’; b, b’;¢; x,y] = A,(c)P(a, b; c; x) O(a’, b’; c; y) 
Oa; b, b’;c;2,y] = as c)P(a, ee pikes 
Oa; b, b’;c;2,y] = V,(a)®®a, a; bb’; ¢; x,y] 
Oa; b, b';c;2,y] = — \Da; b, b’;c,c; 2, y| 
O®[a,b;c,c'; x,y] = V_(b)®®[a; b, b; c, c'; x, y] 
From these transformations twelve interesting expansions of the 
above functions will be obtained by the use of certain theorems 
stated as lemmas in §7 below. 





6. The fundamental transformations for the abnormal func- 

tions 

When we seek a basic analogue for the ordinary hypergeometric 
function F(a,b;c;x+y) an abnormality appears, which of necessity 
introduces another type of function. I term this abnormal. The 
(a+-y)” of the ordinary binomial expansion finds its parallel in basic 
algebra in a product 

[e+Y]n = (e@+y)e+yq)...(e+yq"). 

This is the basic binomial for positive integer n. Its well-known 


expansion is ! 


(n)! N—To Ty hr(r— 
> (r)!(n—r)!" yer, 


introducing the quadratic solitary factor, so that we have the 
abnormal function 
ns b) 
O(a, b;c;|a-+2 > > (a 1) m+ n9) m+n _ gemyngin(n—1) 
( [x+y]) > 2, Gm)! (nyt y"q 


m=0n !(n)! C)m-+n 


0 
S (2), Hb) rtay)(: n— 


From these by means of the lemmas in §7 below we obtain 
O(a, b;c;[x+-y]) = V_(b)O®[a; b, b; c; 2, y;q] 
O(a, b;c;[x+y]) = A,(c)®®[a, b; c,c; x,y; q] . (14) 
O(a, b;¢;[a+y]) = V,(b)A,(c)B®[a; b, b; c,c; 2, y;q] 
The presence of the element g in the functions on the right denotes 
that a quadratic solitary factor g'""-» is associated with y” in the 
general term, and, of course, in every term of the infinite series. 
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7. Lemmas 

Since (a), = (l—q*)(1—gq***)...(l—q**") = (a)(a+1)...(a+r—1), 
so that throughout the paper the terms of factorials are always basic 
elements—thus (a+r) = 1—q*+’—we shall have, in the case of an 
operative factorial (—8@),, 

(—@),2™ = (—m)(—m-+1)...(—m+r—1)a™ 
= (—m),2", (15) 

which vanishes identically for all positive integral m < r. 

Moreover, when m = r, we can, by a transformation, have also 


(—m), = (—Y(r)!q-**, (16) 
and so, by introducing the operator ¢, we obtain 
(—@),(—¢),2™y”™ — (—m),(—n),2™y", (17) 


which vanishes identically for min(m,n) < r. This and the following 
three operations are of the first importance; the usefulness of the 
A, operative functions entirely depends upon them: 


my griatm+n)—inr—1), (18) 
(a+m+n—r), 


{1/(0+a),(d+a),}amy" = {(a+m),(a+n),}-lamy”, (19) 
{(6+-c)}@(a, b;c;x) = (c)-'@(a, b;c+r; 2x). (20) 


{1/(1l—a—0—¢),}a™ wee = (— \ ee 


We have, in addition to 


_ Gh)G(h+0+4¢) _ YS (—9)(— 1044 
Valh) = G Gh OG) ~ a ie Ir grh+0+4), (21) 


its inverse 


(22) 


(h)»(—9),(—¢), inr—Dirk+0+#) (23 
Be (h+),_ = 





(k—h),(—®),(—¢4), gi- —1)4+1(k+0+¢) 


(r)! (k+-r—1),(k+-9),(k+-4),(h), 


s k),(— 6),(— -9) (—),_ r 24 
(r)!(h),(1—k—0—9), = * Ge 
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Further, 

®, 

. | (a), (6) 
G(a)Q a+6+¢) (b—h), —¢). a-241 
G(a+6)G(a+¢) 0)" Bah." rae (b) @ . 

The basic binomial identity is 
If (a-+b)y = (a+b)(a-+bq)(a-+bq?)...(a+-bq), (26) 
(a+5)y mas a’b§ }—s(s—1) 27 
ats ay = > wel — 


It must be noted that the order r, s is important. A change to s, r 
would involve a change in the solitary index from s to r. 
The basic analogue to Vandermonde’s identity is 


If [aly = (1—g*)(1—q**)...(1—q*t-1), _ etc., 


then [a+ bly = > ( a 1b. gq. 


(YN)! (r)! (8)! 


Here also the order of 7, s must be observed.* 


* q*tt- saiaad 


?+s= 


8. The twelve normal expansions 
The detailed analysis for one function will be sufficient to show 
how all the twelve normal expansions may be obtained. 
From (11) we have 
O(a, b;c;x) O(a’, b’;c;y) = V,(c)®®[a, a’; b, b’; ¢; x,y]. (30) 


Replace V,(c) by its equivalent series (21), namely 


«o 
pia r gric+8+9), 


By operating on the general term of ®® we obtain 
(—m)(—)y _ 9e4many(®)m(@’)n(O)m(5’) * 
c+m+n m “wm AN Man | 31) 
HO, % (mM)! Omen ; 
Now (—m), or (—n), vanishes identically for min(m,n) < r. 
We therefore (for such value of m or n = r) replace all m, n by r; 


* Heine, Kugelfunctionen (1878), 99; W. N. Bailey, Generalized Hyper- 
geometric Series (Cambridge, 1935), 69; G. N. Watson, J. of London Math. Soc. 
4 (1929), 4-9; F. H. Jackson, Quart. J. of Math. (Oxford), 11 (1940), 16 (V). 
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then making use of transformation (3) we must form the proper 
solitary factor by putting m = r, n = r in the expression 


qi—2m +r—1)+4r(—2n+r—-1+1(c+m+4n) 
(a),(a’),(b),(6’), 
(r)! (c),(C)oy 


and consequently an expansion 


So we have ger-lyy’, 


O(a, b; c; x)M(a’, b’;c;y) = » eaunren Mr ryt gre) 


x @®[a+r, a’+r;b+r, b’+r;c+2r; x, y)}. (32) 


As similar methods only are needed in order to obtain other 
expansions, it is necessary merely to state the expansions, and in 
doing so, I follow the order of those in Burchnall and Chaundy’s 
paper.* Each expansion is followed by its inverse. 


Dla; 6, b'5¢,¢73 x,y] = 2 a rehet atygie#) x 


O(a+r,b+r;c+r;2)O(a+r,b’+r;c’+r;y), (33) 


(a),(6),(6'), 
(r)! (c),(c’), 


XO%a+r;b+r,b'+r5c+r,c'+r;a,y], (34) 


(a, b; c; x) (a, b’;c’;y) = > (—y> ary’g@ tire) x 
r=0 





O®a, a’; b, b’;c; ay] = Z (— per ee " atytgretine-) x 


Te (35) 


O(a, b; c; x) @( (a’ i, 3c; y) = ph tee Dr gryrgrte+t-D > 


spf +r;b+r,b’+r;c+2r;2,y], 


PY! a; b,b';c;2,y] = 2: ay ee a arytgretinrD x 


xe, bnekbie+i eee (37) 


* J. L. Burchnall and T. W. Chaundy, Quart. J. of Math. (Oxford), 11 
(1940), 253 et seq. 
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-(b),(6’),(c—a@), r% yar+trir—1) 
NO, % * 


O%Ya+r;b+r,b’+r;c+2r;2,y], (38) 


M(a,b;c;x)D(a, b’;c;y) = pi (—_y @ 
r=0 


Oa; b, b’;c;2,y] = - oa Mr vy get) x 
r=0 ; - 


xX O®fa+r,a+r;b+r,b’+r;c+ 2r; x,y], 


: (a),(6),.(6’) ae 
O®)a,a;b,b’;c;x2,y] = (—) tyr ter-e x 
/] > nites 
x OYa+r;b+r, b’+r; c+ 2r; x,y], 
’ ~~ (a).,(b),(b’) ieee 
MW . e, "O° 2, sie —_ ms 2r r r Ea) ragret+irr Ny 
[a;b,b’;c;2,y] 2, | ) yi e+r—ihiele yd 
ere b’+r;e+2r,c+2r; x, y], 


O° a; : > (4) ie) r ry rytgretine— 1). 
viata y|, 


Oa, -¢ c’ at a 5 nes 


r=0 
netet tts b+r;e+r,c’+r; a, y], 
O° la; .b:¢.¢ ‘'32,y ee “Ve Jor (b),. vege) x 
rac 4 Teel 


x ay +2r,b+r;e+r,c’+r;a,y]. (44) 


This completes the expansions in which all the series are of normal 


form. 


9. The six expansions for abnormal series 

By using the operators V,(b), A,(c), V,(b)A,(c) and their series 
equivalents as given in (21)-(27) we can obtain the following six 
expansions. Each expansion is followed by its inverse. 


~ b 
D(a, b;c¢; [w+y]}) _ > moe ary'grot-) x 
r=0 ‘ 2r 


x @O%7a+ 2r;b+r,b+r;c+2r;x,y;q], (45) 
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©a; b,b;c;2,y;q] = > (—y an atytgrtined) » 
0 er 


x @(a+ 2r,b+-r;c+ 2r;[x+y]), 


%a,dicslety) = > (Yep tye ern 


ie b+2r;c+2r,c+2r;x,y;q], (47) 


r=0 


Oa, b; C,C; 2, ¥; q| ay i ary" grctie— Dy 
) Cc) 


ox ( 
p< ee (48) 
“ b),(c—b), 
O(a,b;c;[a+y]) = _(@)xr(),( ~aytgretie—-D x 
ety] z (I (c+r—V),(C)y 
x O°[a+ 2r;3b+r,b+r;c+2r,c+2r;x,y;q], (49) 
‘ ~ (4) 2,(6),(c—b) 
O° a;b,b;c,¢;x2,y;q| = —y—— PT _ gtytgro+r-l) x 
[ yg = > (-) mite, 7"! 
x @O(a+ 2r,b+r;c+2r;[a+y]). (50) 
It is to be noted that, where the element q follows x, y in the 
functions ®®, ®®@, etc., the functions are abnormal. A quadratic 


solitary factor appears in each term. In the general term of such 
functions we always have x’y* multiplied by g***-». 


10. Two exceptional expansions 

Before giving these expansions, it is necessary to introduce a basic 
analogue of Taylor’s theorem due to the present writer.* 

If f[x+h] be an absolutely convergent series 


> Alz+h}, 
in which [a+h], = ere, ee 


oO 


then flx+h] = ~ 2 ge- (DV) (51) 


—f 9. we 
a(l—q) 2x(l—q) 





in which D denotes the basic operator 


operator is the basic analogue of d/dx and is such that 
a{D} = (0)(@—1)...(@—r+1). (52) 
* F. H. Jackson, Messenger of Math. 39 (1910), 26-8. 
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If we apply this general theorem to the special basic hyper- 


eometric function 
& 2, (a, b; c;[x+h}) 


and replace the arguments z, h by (i) x+-y, —ay; (ii) e+y—xy, +2y; 
we have the following two expansions: — 
> (2),(b) . 
O(a, b;c;[x+y | —ary]) = (—)r Sher gryrgit—) x 
: = DY oreY2 


x@(a+r,b+r;e+r;a+y), (53) 


O(a, b; ¢;[e+-y—azy | +2y]) = 3 mt - aygie x 
r=0 of 


x@O(a+r,b+r;ce+r;at+y—ay). (54) 
I term these exceptional expansions. They are basic analogues of 
Burchnall’s and Chaundy’s results (44), (45). The unusual notation 
[x+y | —ay] is necessary to avoid confusion with «+y—ay: thus 
[z+y | —2y], 
= (w+y—xy)(x+-y—axyq)(x+-y—axyq’)...(e+y—axyg), 
[e+y—ay | +2y], 
= (e+ y—ay+ay)(e+y—ary+xyq)...(e+y—axy+axyq'*), 
showing that [x+y—ay | +2y] is not (x+y) except when |q| = 1. 
Moreover, in 
M(a+r,b+r;ce+r;a+y) and O(a+r,b+r;c+r;24+y—2zy), 
we have ordinary basic hypergeometric series in which the variable 
x+y appears in the general term as an ordinary binomial (x+-y)” 
and «+y—zy appears as a trinomial (x+y—zy)". Such series as 
these, in which we have factorials of basic form together with ordinary 
binomial forms, make it impossible to obtain basic analogues of 
Burchnall’s and Chaundy’s expansions (46)—(51). I therefore give the 
name exceptional to the series to avoid confusion with abnormal series. 
I propose to consider generalization of Burchnall’s and Chaundy’s 
expansions (52)—(55) in a second paper dealing also with the basic 
confluent functions. 


11. Duplication formula 
The following special case of the expansion (37) seems interesting 
in connexion with Vandermonde forms. 


If in Oa; b,b';c; x,y] 
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we replace y by xq” and apply the basic Vandermonde theorem, ®® 
reduces to a basic hypergeometric series O(a,b+-b’;c;x), and the 
theorem (37) becomes 


, ~* (c—a),(a),(6),(6'), 
O(a,b+b';¢;2 ames r girgratd)+ine—D) x 
= > leroy? 
x O(a+r,b+r;c+ 2r;x)O(a+r, b’+r,c+2r,aq’), (55) 
which is a quasi-addition formula for a pair of numerator elements. 
In the same way we can obtain an inverse formula using the 
expansion (38): 


, x b04 b’ )y (a a), 
O(a, b; c; x) O(a, b’;c; xq’) = ne yr (C—4), (6),(6"),(@), Pr gp 2rgr(a+b)+inr—l) y¢ 
yo(a, W'seiea") = > (Yate 
x @O(a+r,b+6'+2r;c+2r,x). (56) 

From this, by giving special values to the parameters, we obtain 
q-generalizations of certain theorems first given by Bailey.* 

In (56) put c = 6+-b’; then we have 

(a,b; b+6'; x)@(a, b’; b+-b’; xq”) 
=> (4), (0), (6° )r (6+-b'— @)r p2rgria+b)+4rr—) »¢ 
a (r)! (bb ‘)(b+-8' Jor 


r 





(a+r) , (a+r) 
x {1 (I) Cn t (57) 


which, by well-known expansions and properties of factorials, gives us 


O(a, b; b+b'; x) O(a, b’; b+-b’; xq’) 


rat) Si @xlO) (0 (6-+0'—A), oe 
= (2-0 De) GED ELD )gg TAMA L ONE, (68) 


These reduce to Bailey’s theorem in the limiting case |q| = 1. 


12. The q°, g? partial differential equations 
The basic analogues of p, q, r, 8, ¢ which are generally used to 
represent 2022 eae @ 
aa” By’ Bat” Baby? dy? 
are represented by P, Q, R, S, T respectively, where 
. = 1—q? 
" ai .~ aw R= PF, S= PQ, T=. 


* W.N. Bailey, Proc. London Math. Soc. (2) 38 (1935), 377-84. 
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For our function O® we have 
[(a)(b) +-y(b)q**9Q-+ {9° (a+ 1)a+-q%(b)x—(c)}P+ 
a c= 0, 
[(a)(b’) +-2(b’ )g**+4 P+-{q" (a+ 1)y+-9%(b' )y—(c)}Q+ 
+ farytth' +1 gg S + yfgate' +941 9 TO) = 0. 


In factorial forms the above two ery become 


[y(6-+¢+4)(p+b)—(¢)(8+4¢+c'—1)]O® = 0, 
Sa io = 0. 


The equations for ®® are rather more simple. They are 
a y(b)q2*"Q+{q x+q?(a)a—(c)}P+ 
[(a)(b)+-y(b)q*?Q+ {g(b+1)x+-q"(a)x—(c)}P+ 
parygt1 8+ x(agttot1+d — 9°) RID® = 0, 
[(a)(b’) +-a(b' qt? P+ {g4(b' + 1)y +9" (a)y—(c)}Q+ 
+aygt’ +19 +y(yqtt’ +148 9°) T D2) = 0, 


The equations for 0 are 
[ (a)(b) +-{q%(b+ 1)a+q?(a)a—(c)} = +b+1y7—9*t R}IO® — 0, 
[(a’)(b') + fa" + (6+ Dy +9" (a’)y—(c)}Q— 

—agS—yf{gt +g} T]O® = 0. 
The equations for ® are 

[(a)(b) + {q"(b-+1)x-+-g°(a)x—(c)}P+ 
+ {g(b+ ly +e(a)y}q? Qt qe ay ig’? +-gS+ 
+y 2qa+b+8 HT + afagtto+$+1_ gel RIO =Q 


and a second equation formed by interchanging 0, ¢ and 2, y and 
by replacing c by c’. I need not write down the equation. The 
equations for ®® are the most simple. All the equations reduce in 
the limiting case |q| = 1 to the partial differential equations for 
Appell’s functions as given in Bailey’s Tract on generalized hyper- 


geometric series. I do not proposé to discuss these equations in this 
paper, though it is obvious many results of interest may be deduced 


from them. 








ON BASIC DOUBLE HYPERGEOMETRIC FUNCTIONS 81 
13. A basic integral for 0 
I use the notation S d(gz) to denote the operation which reverses 
3 
* a(1—q)’ 
In former papers* on q-definite integrals I showed that 


Po = 1,(6)U(e—6) 
S wh"(1—qu),-4-a(1— gut) _q Agu) = T(c) O(a, bc; 2), 


1 
SS wA(—qu)g_y-10°(1—qu),-¢-a(1— quer) 4 A qu)d(qr) 
0 0 
L(eP,(e— 


eM g(O)T'(d—6) 3,(a, 6; c;d, e; 2) 
P,(e)T,(d) 
in which [, represents the basic gamma function. 
By antinis similar to those in the papers to which I refer we can 
show that 
Lay, ( 
l(c) 


qd 


C—4) para; b;b’;c;2,y] 


1 . 
= S$ w(1—qu),_g4(1—qreu)_,(1— gu) _y d(qu’). 
Su 


It is not possible to find similar basic integrals for the other functions 
0, O®, O®, since trinomials appear and in basic notation there are 
no simple analogues for multinomial theorems. In the above formulae 
(1—g"A)_, =14 0), PM L— 9") 04 

(1) (2)! 

supposed absolutely convergent. 
Also the basic gamma function has the property 

L(n+1) = (l—q™)0,(n). 
For a full discussion of this function I refer again to the papers by 
Littlewood and the present writer to which reference was given 
earlier in this paper. In conclusion, it is obvious that the subject 
offers great field for development in the case of functions with more 
parameters, and in discussion of confluent cases analogous to those 
given by Burchnall and Chaundy. I leave this, however, to what 


I hope may be a second paper. 


wees 


* F, H. Jackson, Quart. J. of Math. 4 (1910), 193-203; 199; American J. 
of Math. 32 (1910), 305-14. 
3695.13 G 
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Conditions for Convergence 

In order that the expansions may be absolutely convergent in the 
threefold multiple series which appear, conditions must be stated. 
I avail myself of the results given in great detail by Burchnall and 
Chaundy* for their series which are limiting cases when |q| = 1 of 
those of this paper. It is obvious that the same conditions will cover 
the series discussed in this paper not only in the limiting case |q| = 1 
but under the condition |g| <1, which condition I impose. The 
presence, moreover, of the quadratic solitary factor will ensure very 
rapid convergence. 

* Loc. cit. 263-9. 








A THEOREM ON INDEPENDENCE RELATIONS 
By R. RADO (Sheffield) 
[Received 12 November 1941] 


1. P. Hall* considered the following problem: Given a finite number 
of abstract sets A,, Ag,..., A,, under what circumstances is it possible 
to select one element a, from each set A, in such a way that a,, ay,..., Uy, 
are different from each other? He found that the following conditions 
which are obviously necessary are also sufficient for the possibility 
of this choice: However one selects any number k of the sets A,, 
say A, A,,,..., A,,, where lok<on; loiyu<vy<... << <0, the 
sets thus selected contain between them at least k different elements. 
The purpose of this note is to investigate what additional restric- 
tions apart from their being distinct may be imposed upon the 
‘representatives’ a, if the obviously necessary conditions for their 
existence are to be at the same time sufficient conditions. In 
particular, it will be shown that in the case when the elements of the 
sets A, are vectors in some euclidean space, we may postulate that 
the selected elements a, are linearly independent. This leads to 


THEOREM 1. Suppose that A,, Ag,..., A, are sets whose elements are 
vectors in some euclidean space. Then, if no group of k of the sets A, is 
contained in a sub-space of k—1 dimensions (k = 1, 2,...,n), tt is 
possible to select one vector a, from each set A, in such a way that 
Ay, Ag,..., A, are linearly independent. 


As a second example consider the case when the elements of the 
sets A, are polynomials in variables 1,,..., t,,, with coefficients in some 
field F. We call polynomials a,,..., a, independent if no polynomial 

w),,) exists with coefficients in F not all zero which has the 


property that f(a,,...,@),) = 0, identically in ¢,,..., t,- 


THEOREM 2. Suppose that A,,..., A, are sets of polynomials. If 
any k of the sets A, contain between them at least k independent poly- 


* P. Hall, J. of London Math. Soc., 10 (1934), 26-30. 
+ The 2”—1 conditions obtained in this way are independent of each other 
as is shown by the example: 
A, = A, =... = Ay = {1, 2, ...5 K—1}3 Agys = ... = Ag = {1,..., 1} 
in which every one of Hall’s conditions is satisfied except the one referring to 
vy, = 13 ve = 23 203 mm = B. 
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nomials (k = 1, 2,..., n) then it is possible to select one polynomial a, 
from each set A, in such a way that a,,..., a, are independent. 


Theorem 2, and in fact, an analogous more general theorem con- 
cerning arbitrary functions instead of polynomials, follows from 
Theorem 1. For, under suitable regularity assumptions, the functions 
A ii(tys-++stin)s Foltas--+s tin)s-++> F(tas---> tm) are independent if and only if 
the k vectors whose components are the elements forming the rows 
of the matrix (2) , are independent. 

ON pd Ky 

2. Let us now pass on to the general case. Suppose that the 
letters a, u, x, y, z as well as a’, etc. denote elements of some abstract 
set S. Aj, Ag,..., By, By,... denote sub-sets of S. Let J be a function 


which associates with every ordered finite system 2, %p,..., %, of 
elements of S a number J(2,,...,,,) which is either zero or unity. 


For convenience we admit the value m = 0 in which case E(e4,..+. Sad 
is written as /(—), and we suppose throughout that J(—) = 1. 

I is called an independence relation, I(2,...,%,,) = 1 corresponding 
to the assertion that the elements 2,,..., 2,, are independent, if 

(i) Lis decreasing, i.e. 

Pete iitgd Elias +< Hes Maen) 

for every m > 0 and every 2..., % (By removing any elements 
from a set of independent elements another set of independent 


m+l1" 


elements is obtained.) 
(ii) Lis commutative, i.e. 


whenever A,,..., A,, is a rearrangement of 1.,..., m. (Independence of 
elements is not affected by changing the order in which they follow 
each other.) 
(iii) J is non-reflexive, i.e. 
Iz, z) = © 


for every x. (Two independent elements are always different from 
each other.) 
(iv) J is distributive, i.e. 


T(Xy,.--5 Um )L(Y05Ymsa) SY L(%p---s ms Yp) 


for every m > 0 and every 2},..., Ups Yys-> Ymar- (If m elements 


m 
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X1,+++) %m are independent and every one of m-+-1 elements ¥j,..., Ym41 
depends on the x,, then y,,..., Y,,4; are not independent.) 

A system of sets A,,..., A, is said to have the property H (Hall’s 
condition) with respect to J if, given any numbers k, »,..., vy 
satisfying 

r< 2 <= &, liy<wy<..<y KN, 


the sum-set A, +A,,+...+A,, contains elements 2,,..., x, for which 


Sets A,,..., A, are said to have the property R (existence of 
independent representatives) with respect to J if an element a, can 
be chosen from every set A, in such a way that I(a,,...,a,) = 1. 

The function J = J, defined by putting 


El G_5002) Sq) = 1 


om 
if and only if 2,..., 2,, are distinct, has properties (i), (ii), (iii), (iv). 
Hall’s theorem, quoted at the beginning, asserts that H and R are 
equivalent with respect to J,. This is a special case of 


THEOREM 3. With respect to every independence relation I, the 
properties H and R are equivalent. 


THEOREM 4. If, with respect to a non-reflexive function I, the 
properties H and R are equivalent, then I is an independence relation. 


The author has not been able to decide which property should be 
substituted for non-reflexivity in order to make the correspondence 
between Theorem 3 and Theorem 4 a perfect one. It is obvious that 
there are functions J, decreasing, commutative, and distributive, 
but not non-reflexive, with respect to which H and R are equiva- 
lent, e.g. the trivial function defined by putting 

Zz) =I 
for every m > 0 and every choice of elements 2,,..., %,,. On the other 
hand, an example will be given at the end of this note of a decreasing, 
commutative, and distributive function with respect to which H and 
R are not equivalent. 

Familiar facts concerning linear relations between vectors show 
that Theorem 1 is a special case of Theorem 3, and it has already been 
pointed out how Theorem 2 follows from Theorem 1. It does not 
seem possible to establish the foregoing results by an adaptation of 
Hall’s proof of his special case, nor is there any connexion between 
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the present note and an earlier paper* which generalized Hall’s 
theorem in a different direction. 

3. We start by proving the almost trivial Theorem 4. Suppose 
that J is non-reflexive, and that H and R are equivalent. Then J is 
decreasing. For let m >'0; I(2,...,%,4;) = 1. Then the sets 

A, = {x,} (l<p<m+l) 
have the property R and therefore H. From the definition of H it is 
clear that then the system A,,..., A,, satisfies H, and hence R. This 
last statement means that 
I(24,...5%_) = 1. 
Thus, bearing in mind that J takes only the values 0 and 1, we find 


that ’ . i. , eae 
Fhe cncei Men) oF Bigs. a 


’ 
m-+1° 


Next we show that J is commutative. Suppose that 
ne Ln) = 1 


for some m > 0 and some elements 2,..., ,,. Then the system of 


for any m > 0 and any 2},..., 2 


sets A,,..., A,, where 4 fix 
a a BS 
satisfies R and therefore H. From the definition of H it follows that, 


for any permutation 
a aan 
ye a 


the system A), A),,..., Ay, satisfies H, and therefore R. This is the 
same as saying that 


ey 
= {2 


I(x geeege ry,,) = 


, 


Thus, for any m > 0 and any elements 2},..., Xp; 
E(B) ,:<-y Bh) = TG j,---5 Zqy)- 
Finally, J is distributive. For let m > 0, 
I(xy,..-, Xm) or I(y, Aig Ym+1) =1. (1) 


Put A, = {a,} (l<n<m), 


Ams1 ae {1 peate Ym} 
Then A,,..., A,, satisfy R and therefore H. Now, since A,,,, con- 
tains m+ 1 ‘independent’ elements, a moment’s consideration shows 
that A,,..., A,,, A»4+ satisfy H. Hence these sets satisfy R, i.e. for 


m? m+ 


* R. Rado, Proc. London Math. Soc. 44 (1938), 61-91. 
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at least one value of jug (1 < pp < mM+1) we have 1(a,..., Xp, Y,,) = 1. 
Then onda 


Dd L(Sp.--2 Lm Yp) S U(X y,---1 Lm) L(Yr---> Ysa) (2) 


p=1 
In view of the fact that J takes only the values 0 and 1, we see that 
(2) holds for any 2},..., Zp» Y1---> Ym41» respective of whether (1) is 
true or not true. 


4. We now proceed to prove the principal Theorem 3. Consider 
a decreasing, commutative, non-reflexive, and distributive function 
I with respect to which H and R are to be taken. The implication 
R > H is trivial. For, if A,,..., A, satisfy R, then every A, contains 
an element a, which is such that I(a,,...,a,) =1. Ifnowl<k <n; 
liyu<wy<..<¥<n, then the set A,,+Ay,+...+A4,, con- 
tains the elements @,,, @,,,..-, @,,, and we have 


2 VE? 
I(a 


yyereo Ay) DS I(a,,...,4,) = 1. 


Therefore A,,..., A, satisfy H. 

On the other hand, let us suppose that the sets A,,..., A, have the 
property H. We have to prove that they satisfy R. If n = 1, then 
H and R are identical in meaning. We may therefore assume that 
n > 1, and that the equivalence of H and R# has already been 
established for any system of n—1 sets Aj,..., Aj,,. Then there are 
elements x, in A, (1 < v < n) satisfying 

IE ae5,.-, Mga) = 1. 


For, since A,,..., A, satisfy H, the same is true for A,,..., A,_,, and 
then the existence of the x, follows from our induction hypothesis. 


Furthermore, since A,,..., A, satisfy H, we can find elements 
Yis-++> Y, in the set A,+A,+...+A, for which 
I(¥,,---5 9%.) = 1. 


The distributive inequality 


1 a= F(a,,..., Sq—a)TGg-:-1 Fa) SD, 1la,,---1 qa Hed 


v=1 
shows that for some index vp (1 < vy < n) we have 
} | Tn-p Yy,) = : 
Case I. Suppose that y,, is contained in A,,. 
Then 2},..., Z,-1, Yy, are in Aj,..., An, A, respectively, and this 
last system of sets has the property R. 
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Case II. Suppose that y, belongs to A), where 1 < Ay < n. 
Putting kwd,, 2B «ih; 
B= Svgn;v 41,A)), 


" 

> — , ~ a . a’ — 

i ©.» —F ? 7. Yy> 
2, =e <v<n;v Aj), 


we see that +s 
. z, isin B, 


2; is in B,, 
I (a;, 2 
Furthermore, B,,..., B, satisfy H. In particular, B,,..., B, satisfy 
H, and hence, by induction, R. There are elements u, in B, 
(1 <v< _n) for which 
E(ttq,.-.,%,) = 1. 


Define r(ug,..., u,,.) as being the number of indices v with 1 < v < 


for which u, = z,. Then 

0 < 7r(Us,...,U,) K N—2. 
Choose %p,...,u,, in such a way that, for the 2,...,z,_;,2, under 
consideration, 7(wp,...,¥,) is maximal. We have 


Fe ee ee eee es 
(3) 


n—1 
I (Ug,..+5 Uns 2%) + > L(tUg,..-5 Uns Zy) 
v=1 


Therefore at least one term on the right-hand side of (3) is positive. 
Case II’. Suppose that I(ug,..., u,,, 21) = 1. 
Then 2}, to,..., u,, are respectively in B,, B,,..., B,, and this last 
system has the property R. Hence A,,..., A, satisfies R. 
Case II". Suppose that J(ug,..., u,,, 2) = 1. 
Then precisely the same conclusion follows. 
Case II”. Suppose that, for some index vy, (1 < », < n), 


E(the,..., Un, Zy,) = 1. 


nm? Vy 


Then I(t, Zy,) 2 I(teg,...) Was Zy,)s 


Uu 7 % 


V1 v;? 
I (tg,..., By. -39 2y,> 
The elements a 


Mm-b “vy? 


are respectively in ,,..., B,_,, B,,, 
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and we have 
1(Ug).20) Uy.) Zy,9 Uy tzrreer Uy) = 1(Ug,---y Uyy--ey Uy) +L, 
which contradicts the fact that r(ug,...,u,,.) has its largest possible 
value. This concludes our proof. 

The example, referred to in §2 above, of a function J which is 
decreasing, commutative, and distributive, but with respect to 
which H and R are not equivalent is as follows. Put S = {1, 2, 3}. 
Define J by putting 

I(—) = I(x) = I(x, y) = I(x, y,z) =1 
whenever the set {x,y,z} contains exactly two different elements, 
and I(a,2%,...,% »,) = 0 in all other cases. It is easily verified that 


oom 


I has the required properties. If we put 


A, = {1}, A, = {2}, A, = {3}, 
then the system A,, A,, Ag satisfies, with respect to J, Hall’s condi- 
tion H. On the other hand, there are no independent representatives 
of A,, Ag, Ag, ie. R does not hold. J is, of course, not non-reflexive. 








DIFFERENTIAL EQUATIONS ASSOCIATED WITH 
HYPERGEOMETRIC FUNCTIONS 


By J. L. BURCHNALL (Durham) 
[Received 11 March 1942] 


1. Introduction . 
In the course of a discussion with Mr. Chaundy regarding the ‘three- 
term’ differential equation 

{ f(8)—a"g(5)+27"h(5)}y = 0 (5 = ad/dz) (1) 
it occurred to me to form the linear ordinary differential equations 
satisfied by certain functions F(px,qx), where p, q are parameters, 
x is the variable of differentiation, and F is a hypergeometric func- 
tion of two variables of a type specified in each case. Except in the 
final illustration, where I extend my method to Whittaker’s func- 
tions, the equations obtained are of the type (1), and it seems worth 
while to record these instances in which an equation, so intractable 
in general, may be solved in terms of known functions. One point 
of special interest, discussed in detail later, may be mentioned here. 
If we take F to be a product of Bessel’s functions of different orders 
and arguments, my method gives an equation of the fifth order 
instead of the anticipated equation of the fourth order. The former, 
though possessing an irrelevant solution, is in some respects a simpler 
equation than the latter, and the irrelevant solution has a simple and 
interesting form. I show in §9 that a similar state of affairs exists 
when we consider products of functions which satisfy a large and 
important class of differential equations of the second order. 

My method is based on the elementary observation that 


. é d 
} P(px,qz)+a5 F (px, qx) = x Ppa, ga) 


o 
) — — 
op x 


or, setting Pz = 6, 17, = ¢, zo = 6, that 


0+ = 8, (2) 
for operations on the functions considered. In practically all the cases 
dealt with the function F satisfies two ‘two-term’ partial differential 
equations 


fi(9, 6) F _— pxg,(8, ¢)F, f2(9, $)F a gxg,(8, ¢)F, 
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where f,, g, are polynomial in their arguments with coefficients in- 
dependent of p, g, x. By suitable manipulation it is possible to obtain 
from these equations an equation in which 0, ¢ occur only in the 
association 6+¢. On replacing 6+¢ by 8 we have an ordinary 
differential equation satisfied by F regarded as a function of z. 

The connexion between systems of partial differential equations 
and ordinary equations may of course be explored in either direction, 
and I am indebted to a referee for drawing my attention to a series 
of papers by Horn, of which (4) is the latest, written from the 
opposite point of view to that taken here. Horn’s aims and methods 
are markedly unlike those of the present paper and the points of 
contact between our work would seem to be slight. 

We may illustrate the method and provide a starting-point for 
further development by considering the function 


z= (1—pax)*(1—qa) = ,F(b; px) ,F,(0,, qx), 


where the notation is that usual for hypergeometric functions. Since 
(1—zx)~® is a solution of the equation 


dy = x(5+5)y, 


it is evident that z satisfies the partial differential equations 


0z = px(6+b)z, oz = qu(d+b’)z. 
From these, by addition and use of (2) and (3), 
z = (0-+¢)z = px(S-+b—9)z-+qx(8-+b’—0)z 
= af p(6+6)+9(5+0’) jz—pqu*[ +b’ + 6+]. 
Thus z satisfies the ordinary equation 
Sz—a[ p(5+b)+4(8-+0') e+ pgx*(5-+b+b')z = 0. 
Suppose now that the equations (3) are replaced by 


Of(6+-¢)z = px(6+-b)g(6+-¢)z 
of(0+-¢)z = gu(d+b')g(0+-¢)z 


Then, by a similar manipulation, 


5f(5)z = pa(d+b—¢)g(8)z+-g2(5+-b’ —8)g(8)z, 


8f(8)2—[px(8-+b)-+gx(3-+b')}e-+ pxgg(8)z-+qx0g(3)z = 0. 
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Operate now on (6) with f(6—1) and observe that 
S(8—1)pxdg(5)z = px f(d)dg()z 
= pxg(d)qx(p+b')g(d)z 
= pqx*(d-+-b')g(8)9(8+ 1)z, 
with a similar result for the term in g. Thus, performing the operation 
stated on (6), we see that 
a function z which satisfies the partial equations (5) also satisfies the 
ordinary equation 
8 f(3) f(S—1)z—a[ p(5-+b)+-9(8+6') | f(8)g(8)z+ 
+ pqu?(3+-b+b')g(8)g(5+1)z = 0. (7) 
There is of course nothing in the argument which precludes the 
existence of solutions of (7) which are not solutions of (5). 


2. The equations of F[a;b,b’;c; px, qa] 


The function 
- (a . +- n (b ae (6' ) vt ( px ym (qa)” 


; aa 
Fa; b, b’;c; px, qa] = » static} 
ot adi m+n 


m=0 n=0 
satisfies the partial differential equations 
0(6+-¢+¢e—1)z = pa(0+b)(6+¢+4)z, 
$(0-+¢-+e—l)z = gx(p-+b')(0-+4-+a)z, 
which are of the form (5) with 


f = 6+¢+c—1, g = 0+¢+4. 


Thus, applying (7), 


Fla; b,b’;c; px, qa] is a solution of the differential equation 
§(8-+-c—1)(8-+c—2)z—a[ p(8+b)-+-¢(8+0’) |(6+-e—1)(6+a)z+ 
+-pgx(3-+b-+b')(8+-a)(S+-a+1)z = 0. (9) 

If in this we make the substitution z = x!-‘z’ and suppress’ the 

accent, we obtain the equation 

8(8—1)(8+ 1—c)z—a[ p(+b+ 1—c)+9(5+b’+ 1—c) ]8(8+a+1 
+pqu?(d+b+b’+1—c)(8+a+1—c)(8+a+2—c)z = 0. 

The singular points of this equation are 


2 = 0, 0, pp, ¢“. 
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The substitutions 


— wyc-a-1 — 7-1 — m-l — q-l 
ay = OO, Hy = ZI", = PH =F 


transform it into an equation identical in form with parameters 
a, = 6+b’+1—c, b, = b, bj = b’, ce, = b+ b’+1—a, 
while the change of variable x, = p-!—z also preserves the form of 
(10) with new parameters defined by 
a, = 1—b, b, = 1—a, b, = b’, c, = c+ 1—a—b, 


Po = Ps Wo = PYG—P)*: 
Thus all the singularities have the same character, their exponents at 


0, 00, i q* 
being 0, a—c+l, 0 0 
1, a—c+2, 1 1 

c—l, b6+06’+1—c, c—a—b, c—a—b’. (11) 





We note that, in accordance with Fuchs’s relation, the sum of the 
exponents is 6. 

Since in (10) the operator 6 is a factor of the middle as well as of 
the leading term, the unit difference of the exponents 0, 1 at the 
origin will not give rise to a logarithmic solution of the equation. 
There will indeed be two series solutions led respectively by a con- 
stant term and a multiple of 2, or, alternatively, two independent 
solutions led by a constant term. Since all singularities of (10) have 
the same character, the equation will have no logarithmic solutions. 
This, evidently, is equally true of (9). 

Before investigating the solutions of (9) valid near the origin we 
may remark that, if in (10) we set 
p=l1,q=y", 6+0'+1—c=a, a+1—c= 8, b’ = f’, 2+-b’—c=y, 
we obtain the equation given by Appell* as that satisfied by 
. Fa; B8,B’;y;2,y] when regarded as a function of 2 only. The 
relevance of this result to the first formula of (13) below will be 


obvious. 


3. Solutions valid near the origin 

The table of functions F® given by Appell} furnishes sixty solu- 
tions of (9). Between any four there is a linear relation with coeffi- 
cients independent of p, g, x, and there are, moreover, triads between 


* (8), 76. 
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which no such relation exists. There is thus a presumption that (9) 
can be completely solved in terms of functions F®, and this is in 
general the case. We may indeed select from the table as solutions 
valid near the origin the following functions: 


2, = Fa; b, b’;c; px, qx) 


v 


= > @n(b )n (qx)"F(—n, b; 1—n—b’; p/q) 
—~ 1! (c), 


oO 


= > ()n(O-+6'), (qx)"F(—n, b;b-+b'; (q—p)/q), 


n! (C)n 
where F is the ordinary hypergeometric function ,F,, and the coeffi- 
cient of x” is a polynomial in p, q of degree n;* 
214 = v1 F(6+b'+1--c;a+1—c, b’; 2+-b'—c; px, p/q) 
(6+-b’+1—c),(a+1—c),, 
n! (2+6’—c),, 
x F(b+6’+n-+1—c, b’;b’+n+2—c;p/q); (13) 


gi-e (pa)" x 


n=0 


and 


245 = v-°F%(b+-b'+ 1—c; b,a+1—c; 2+-b—c; g/p, gx) 


a 


ae > (b+-b'+ 1—c),(a+1—c), 


aA" 


x F(b+6'+n+1—c,b;b+n+2—c,q/p). (14) 


— 
n=0 


If in (14) we employ 
F(b+6b’+n+1—c,b;b+n+2—c;q/p) 
. T(b+n+2—c)l(e—b’—n—1) 
— b+b’+n+1—c ~ \ i eS eae s 
= (—PID) r—b’)r) 
x F(b+-b’+n-+1—c, b’; 2+-6’—c; p/q)+ 


gpl Ont 2—o)l'6'+n+1—c) 
ce oS anda : 
Pid I'(n+2—c)l(b+6’+n+1—c) 

x F(c—n—1,b6;c—b’—n;p/q), (15) 
we obtain a further solution 


= 0S OH eth te 


“15 n! (2— 
fon. n! (2—c),, 


(qu)" F(c—n—1, b;c—b'—n; p/q). 
(16) 
It is natural to adopt z,, 2,4, 21; a8 appropriate solutions when 


* Cf. (4), 450. 
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|p/q| <1. Solutions appropriate when |p/q| >1 may be obtained 
from these by permuting p with g and 6 with b’. We may add that 
the formulae connecting contiguous hypergeometric functions may 
be employed to verify directly that the power series in x occurring 
in (12)-(14) and (16) are solutions of (9). 


4. Some special cases 

We have seen that (9) may in general be solved in terms of func- 
tions F®, but to this general statement I note two exceptions of 
some interest. 

If p = q, which we may note in passing is a singular line of the 
system (8), we may without loss of generality assume p = g = 1. 
The point « =1 is then a confluent singularity of (9) and (10): 
investigation shows it to be regular, and it remains true that solu- 
tions valid in its neighbourhood are free of logarithms. On inspecting 
the solutions obtained in §3, we find that z,, and z,, are, save for 
a constant factor, identical. The equation (9) may now be written 


[(8+e—2)—ax(8+-a)|[8(8+c—1)—a(5+-a)(8+-b+5’)]z = 0 
with the obvious solutions* 
Z, = F(a,b+6';c; 2), 


zg = x1-¢F(b+-b’+-1—c,a+1—c; 2—c;2), 


which are the degenerate forms of z,, 2,4. 
Since (8+c—2)y = x(5+a)y is satisfied by y = x?-*(1—2z)°-4-2, 
any particular integral of 
8(8-+-c—1)z—2(8+a)(6+6+6’)z = Cx?-“(1—2)*-4-2, 
where C is any convenient constant, provides a third solution of 
(9), and, by variation of parameters, we obtain for this third solution 


x 


2 = [ (1—t)-22,,(t)z9(xx) —z4(a)zg(t)} at. 


. 


Another case of failure is at first sight more difficult to understand. 
We observe that, if a+-1 = c, the solutions z,,, 25, 215 all reduce to 
xi and we require an additional solution of (9). It is slightly more 
convenient to consider the equivalent equation (10), which may now 
be written 
(8—1){8+1—c—a[ p(5+6+ 1—c)+9(8+6’+1—c)]+ 
+ pqa?(d+-b+b’+1—c)}8z = 0. (17) 
* Cf. (1), 23, formula (25). 
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An obvious solution is z = 1: to obtain others write 6z = 2°-!u, when 
(17) becomes 
(5+-c—2){S—a[ p(8+b)+-9(5+0’)]+ pqa?(5+-b+b’)}u = 0. (18) 
By (4) one solution of (18) is u = (1—px)-°(1—qax)~”, leading to 
z= [ xe-?*(1—pax)-°(1—qau) dx 


= (c—1)-1¢¢-1 FO (ce—1; bb’; c; px, ga), (19) 


i.e. back to an appropriately modified z, as a solution of (9). Ele- 
mentary methods give as a second solution of (18) 


u = (1—pa)-*(1—qa)-”’ [ x1-*(1—pa)’(1—qa)"’- dx 


leading to a series solution of (10) which does not appear to be 
hypergeometric. 

To explain fully the reason for the present case of failure we must 
recall that the function F™ satisfies not only (8) but also* the 
further equation 

Az = q0(¢+b’ )z—pd(6-+-b)z = 0, (20) 
which is stated by Appell* to be a consequence of (8). It is upon 
the system of three equations that his discussion of the sixty func- 
tions F® is based, and ke has apparently failed to observe that the 
elimination he performs is impossible when a+-1 = c. We may state 
the matter briefly thus. From (8) we may deduce 

(6+¢+c—2)Az = 0 = (6+¢+a—1)dz. 
Thus either Az = 0 or a+-1=, and in the latter case the equations 
(8) have the solution z = (px)-%b(p/q), where ¢ is an arbitrary func- 
tion. This is evidently not a solution of (20). 

If in (8) we write 

(0+¢+c—l)z = (6+¢+a)z = u, 
then the equations for u take the form (3), and so, sufficiently, 

(8+c—l)z = u = (l—px)-(1—qa)", 


whence, in agreement with (19), 


re id [ xe-*(1— px)-°(1—qa)—' dx. 


Thus, if a+1= cc, the system (8) possesses only two solutions which 


* (1), 53. 
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are independent functions of x, and the necessity for an additional 
solution of (9) or (10) is explained. 


5. Extensions and particular cases 
The hypergeometric function of two variables 


(a, )in+n: (a, r/m+ nl) (6 i* m n 
F (px, qx) “22 m!n! (pxy"(qz)", 


bes +n* (Cy) ont TH 


where, for convergence, s > +, satisfies (5) with 
f(0+¢4) = “TT (0+¢+¢—1), g(0+¢) = [1 6+¢+a), 


and is thus a solution of an ordinary equation of form (7) of order 
2s+-1. A complete discussion of all the solutions would appear to 
require the construction of a table corresponding to Appell’s table 
of Fs, 

In particular the confluent function 


, b) ah 
,(b, b’; c; px, qx) => > aie = -(pa)™(gar)” 


satisfies the equation 
§(8-+-e—1)(8+c—2)z—al p(8+6)+9q(5+0’) |(5+e—1)z+ 
+pqx(8+b+6')z= 90 (21) 


whose other solutions are readily deduced from the results of §3. 
We may, however, note that, if c = b+-0’, then (21) may be written 


(5+-b+b’—2)[8(8+-b-+b’—1)—a{p(5+b)+-9(5+6’)}+-pga?}z = 0, 
and that, as originally pointed out to me by Mr. Chaundy, the 


functions 
,(b, b’; b+-b’; pa, gx), x1, (1—b’, 1—b; 2—b—’;; px, gx) 
are solutions of the second-order equation 


8(8-+b+b'—1)z—x[ p(8+b)+9(8-+6')}e-+pge%s = 0. 


6. The equation of ,/\(c, px) »F,(c, qx) 
The function z= oF (c, px) oF (c, gx) 
satisfies the partial equations 
0(6-+-c—1)z = paz, $(¢+c—1)z = gaz. 
By addition we have 


5(5-+-c—1)ze—(p+q)xz = 26¢z. 


H 
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Now 
(8+-c—1)(8+c—2)20gz 
= 2(8-+e—2)6(0-+$-+e—1)z 
= 24(8-++e—2)(6$--paxr)z 
5-+c—2)(6¢6-+-qu)z 
+c—2)0¢z+ pxd(d+-c—1)z+-gx6(5-+ce—1)z 


20( 
26( 
5(5 
_p (8-+-c—2)-+(p+-q)x]6d2+ 2pgqu*z, 
i.e. [(8-++-e—2)(8+ 2c—2)—(p+q)a]206z = 4pqu*z. (24) 
Combining (23) and (24) we have 
[ (8+ 2¢—2)(8-+-c—2)—(p+q)a][8(8-+ce—1)—(p+q)a]z = 4pqxz, 
i.e. 
8(8-+-c—1)(8+c—2)(6+ 2c—2)z— 

—(p+-9)x(26+ 2c—1)(8+-e—1)z+(p—g)*a*z = 0, (25) 
as the differential equation satisfied by )F,(c, px) ,F(c,gx). Other 
independent solutions are 

xi—¢ )F(c; px) oF, (2—e; ga), xi oF (2—c, pa) oF, (c; gx), 
x22 0 (2—c; px) oF, (2—c; ga). 
If the equations (22) are replaced by 
6(0-+c—l1)z = pxf(0+¢)z, $(d-+-c—1)z = ga f(0+¢)z, 
we obtain in place of (25) 
8(8-+-c—1)(6-+ce—2)(8+ 2e—2)z— 
—(p+q)a(5+-e—1)(28+ 2e¢—1)f(3) 
+(p—q)*x S9(8)f(6-+ Ve = 0. (26) 
In particular, taking f = (@+4¢+a)(6+¢-+5), we have the equations 
satisfied by Appell’s function 
which is therefore a solution of 
8(8-+c—1)(6+-c—2)(8+ 2c—2)z— 
—(p+q)a(d+e—1)(26+ 2c—1)(8+-a)(5-+-b)z+ 
+(p—q)*x*(8-+a)(8--a+1)(8+6)(8+b+1)z = 0. (27) 


Other solutions of this equation are* immediately obtainable. 


* (1), 52. 
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Returning to (25), we recall that 
of (c, —fa*) = P(e)(}x)**J,_(2), (28) 
where J is Bessel’s function, and, making the necessary changes of 
variable and writing c—1 = p, we obtain as the differential equation 
satisfied by the products J,,(ax)J,,,(bx) 
5(8—2)(5°—4y:7)z+ 2(a?+-b?)a*8(5+ 1)z+ (a2—b?)2axtz = 0. (29) 
If in this a = b = 1, then it may be written 
(3—2)[3(3*@—4y2)z-+-42°(8-+1)2] = 0, 
and we recognize the inner bracket as the differential equation 
satisfied* by Ji(2). 
7. The equation of ,F,(c; px) F,(c’; qx) 
The partial equations satisfied by this function are 
6(6-+c—1)z = paz, o(p-+c’—1)z = gaz. 
Thus (8+-c+c’—2)6¢z = dpxz+Ogrz = (pxd+qzx6)z, 
and Uz = [6(6+c—1)—(p+q)a]z = (20+c—c’)dz, 
U’'z = [8(6+c¢’—1)—(p+q)a]z = (2¢6+¢’—c)6z. 
6(5+-e—1)(5+-e’—1)z = (8+ce’—1)[06+0(0+c—1)}z 
= (6+c’—1)0dz+ px(d+c’)z, 
$(5+c—1)(8+c¢’—1)z = (6+c—1)0dz+-qau(d+c)z. 
By addition 
Vz = 6(6+c—1)(6+-c¢’—1)z—pa(d+c’)z—gqa(i+c)z 
= (28-+c-+c’—2)6dz. 


Also 


Thus 
(8+c+c’—2)Vz = (26+c+c’—2)(pad+gqaA)z, by (31), 
= px|2¢(¢+ce’ —1)+(26+c—c’)d+ 26]z 
+-ga[ 20(0+c—1)+ (26-+-¢’ —c)0+ 26]z 
= 4pqu*z+pxUz+qxU'z+ 206(5+c+c'’ —2)z, 
by (31) and (32), 


Wz = 8(8+ce—1)(6+c’—1)(6+e+e’—2)z— 
—pa| (8+e’)(8+ce+e’—1)+8(8+e—1)}z— 
—gx| (5+-c)(8+-e+e’—1)+8(3+-c¢’—1) |z+ (p—gq)?a*z 
= 204(5+-c+c’—2)z. (34) 
* (2), 146. 


3453'74B 
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From (33) and (34), 
(28+c+c’—2)Wz = 2(6+e+c’—2)Vz, 
a relation in which 6, ¢ appear only in the association 0+¢ = 6, 
and, on rearrangement, we find that 
z = oF(c; px) F,(c’; gx) ts a solution of 
8(8-+c—1)(6+e¢’ —1)(8+e+ce’ —2)(6+ e+ 4c’ —2)z— 
~2(p+-q)a(S-+$e+4e'(3-+ Jet $e’ —$)(8+-4e-+4e’—Nz+ 
+3(p—q)a(c—e’)(c+-c’ —2)(5+ e+ $e’ —})z+ 
+ (p—q)?x?(5+4c+ }c’+1)z = 0. (35) 
Other solutions are evidently 
~¢ | F,(2—c; px) oFi(c’; 92), al oF (c; pa) of, (2—c’; ga), 
y2-e-c 0B (2—c; px) oF,(2—c’; gz), 
but there is no obvious solution corresponding to the exponent at 
the origin 2—4c—}c’. Reserving this point for discussion until (35) 
has been transformed into the form appropriate to Bessel’s functions, 
I observe that, if for brevity we write (35) as 
fo(5)z—2(p+-9)% fi(8)2+-3(p—Q)x fo(5)2+ (p—q)*x* fa(8) 
then, * a reasoning — to that of § 6, 
So(8)2— q)x f,(8)(8+-a)(5+-b)z+ 3(p—q)x fo(8)(8-+a)(8+-b)z+ 
| ae egg a+1)(8+b)(8 su ~0 (36) 
is satisfied by z= F(a; b;c,c’; px, qx), 
and that three other solutions are immediately obtainable.* 
If in (35) we set c’ = c, it may be written 
(8-+c—1)[8(8-+-c—1)(8+c—2)(5+- 2c—2)z— 
—2(p+q)x(5+e—1)(+¢e—3)z+ (p—q)*x*z] = 0, 
where the inner bracket is recognized as (26). A similar substitu- 
tion in (36) leads to (27). Again in (3¢€ 
factorizes as 
[8+4ce+ dc’—2][8(8-+c—1)(8+-c’—1)(8+e+ce’—2)— 
—4x(3+a)(8-+6)(8-+ fo-+Je’)(8-+4e+-Je’—})z = 0. 


(2¢ 
) 


set p = gq = 1, when it 


* (1), 52. 
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Since the equation possesses only one solution in the form of a power 

series led by a constant term, we evidently have 

a,b, 4c+-4c’, $c-+-4c’—}4; 4x 
c,c’,c+c’—1 

of which (3) 124 (66) and (2) 147 (1) are limiting cases. 


FO(a;b;c,c’;x,2) = I 


| @9 


8. The equation of /J,(ax)J,(bx) 
If in (35) we write 
x=2'* p 1g? —}b?, c—1= yp, c’—l=p, 
(2’)-#-P2’, 
then, recalling (28) and suppressing accents after transformation, we 
find that 

The products of Bessel’s functions J,,,(ax)J,,(bx) are four indepen- 

dent solutions of the differential equation 
= (8—2)[54—2(u?+ p?)d*+ (u?—p?)?]z+ 2(a?+-b?)a*8(8+ 1)(8+ 2)z— 
— 2(a?— b?)(u?— p?)a?(5+ 1)z+ (a?—b?)?24(8+-4)z 
= @, (38) 
If in this we set a = b = 1, the resulting equation may be written 
(5—2)[{54—2(u2+ p?)8*+ (u2—p2)%}z-+ 40°(8-+4+ 1)(8-+2)z] = 0, (39) 
where the inner bracket furnishes the equation of the fourth order* 
satisfied by the product of two Bessel’s functions with the same 
argument. Again, on setting p = p in (38) it may be written 
8[5(5—2)(8°?—4:?)z+ 2(a?+-b*)a8(5+ 1)z+ (a2?—b?)*ar4z] = 0, 
and we recognize the bracket as (29). 

There remains the problem of a fifth solution of (38). We begin 
by comparing that equation with the equation of the fourth order 
whose only solutions are J,,,(ax)J_,(bx). On applying the methods 
of (2) 146 the latter is found to be 
Az 

(p? 4% 2(q pid p2—p?)9} + 2x*(a2-+-b%)(5-+1)(8-+2) e+ 
)ac®[34— 283 — 2(u?+-p?)3?-++ 4(u?+ p?)5+ 3(u?—p?)?]z+ 
)ar*[ 2(a?-+-b?)5(8+ 1) —3(a?—b?)(u?—p?) |z+ (a2?—b?)fx8z 
(40) 
* (2), 146. 
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a result which, in this or an equivalent form, is presumably to be 
found in the literature, though I have failed to discover it. As a 
check upon the calculations we may set in turn p = p and b = a, 
when we recover (29) and (39). 

It will be observed that, whereas Az = 0 is a ‘three-term’ equation 
of type (1), the equation A, z = 0, though free of the irrelevant solu- 
tion, is a ‘four-term’ equation and thus in one respect less elementary 
in character. Since A annihilates all functions annihilated by ‘A,, the 
latter operator is a ‘post-factor’ of A, and inspection of the terms 
of highest order in the operators shows that 


A = {(a?—b?*)a?—(u?—p?)} YS—2—u(z)}A), 


where u(x) is an undetermined function of x, which, by a further 
comparison of coefficients, is shown to be zero. We thus have the 
identity of operators 
[(a*—0®)2*—(2—p2)JA 
= (6—2)A, 
= (5—2)[{(a?—b?)a?— (u?—p?)}84— 2(a2—b?) x88... ]. 
(41) 
Now, for brevity, set J,,,(ax)J,,(bx) = J,, J,,, and let functions 
W.u,xp) be determined to satisfy the equations 
= ie at 0 = 0 
> 5. poUsu,+p = 0 
> 8 J. es Shi, == @ 
[(a®—b2)22—(y2 — p2)e-2 > P(S.p Fup) dUsy,.9 = constant. 


Then z = > a. pwc will be a solution of Az = 0. 


Now, from (41), the 6-Wronskian of the four functions Jn I5, is 
a constant multiple of (a?—b*)a?—(u*—p*), and we find without 
difficulty that the 5-Wronskian of 7“ three _ J, J. » Jon J,, 
J_,J_, is a constant multiple of J_,(bx)dJ_,,(ax)—J_ (axyeJ —p(5x). 
Thus for the solution of (42) we have, Me 


{d= p(bx)dJ-..,(ax)—J;,, (ax)dJ-,(bx)} 
[(a?—b2)a = ‘—2)P 

where the sign of the expression on the right is to be positive if the 

suffix on the left contains like signs, and otherwise negative. 


Buy, = ee ? (43) 








ON HYPERGEOMETRIC FUNCTIONS 
In particular, writing 5, = td/dt, we have 


Usu, ‘e 


aL p(bt) : J. = J_, (at), J_,(bt)| dt 
7 =| —b?)t?—(u?—p*) 
= (b?—a?) PB apeioit (u?—p*)|“[ J_,(bx)5J_, (ax) —J_,(ax)8J_,(bx) | 


* [_,(bt)8 J, (at) —dI_,,(at)83 J_,(bt)] dt 
— fa {ee (2—p) 
(62—a2)"[(a2@—b)x*@—(u2—p®)] > T_,(b22)8I_,(ax) —J_,(ax)8J_,(b2x)} 


x 
+ (b?—a?)- | en 
with similar formulae for the other w’s. 
On multiplying the w,,, ., by J,,,J,, and summing, the sum of the 
terms arising from the integrated parts of w_,,,, is zero and, if we 
suppress the unnecessary factor —* we find that 


Jy, at bt) dt 
z= > +d,,,(ax)J. lt) f eee, 
r | J,(ax), Jj,(at) | | J, (bx), Ps, (bt) | ed 
|J_,(ax), J_,(at)|*|J_,(bx), J_,(bt) 
is the required fifth solution of (38). Since the lower limit ae integra- 
tion is arbitrary, (44) may, moreover, be regarded as the general 
solution. It is an interesting and not too laborious exercise to show 
directly, by differentiation under the sign of integration, that (44) is 
a solution of (38). 


i.e. Sam (44) 


9. An equation of the fifth order associated with two dif- 

ferential equations of the second order 

The results of the last paragraph suggest the following problem. 

Let D = d/dx and 

(D?—I)y = 0, (D?—J)z = 0 
be two linear differential equations of the second order in canonical 
form. Let them have respectively pairs of independent solutions 
Y,(x), Yo(w), and z,(x), z(27) such that 
Y2 Dy, —Y; Dyg = % Dz,—2, Dz, = 1, 

and let 


. 
Y » —_ Y;(x), y;(t) " Z : 
is Y(x), Ys(t) ana 
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We note that 
Fiz) = Z2(z,2)=0 
[DY (a, t))-¢ = [DZ(z,t)|-. =1 

D*Y (x, t) = I(x)Y (a, t), D?Z (x,t) = J(x)Z(a, t) 
Let now WwW = ( Y (a, t)Z(a, t) dt. (46) 
We seek a linear differential equation satisfied by W. Since the 
lower limit of integration is arbitrary, any product of solutions of 
the two second-order equations will be a solution of the desired 
equation. We anticipate therefore an equation of the fifth order. 


Now, bearing in mind (45), we have 


DW = { (Z2DY+YD2) dt, (47) 


(D?—I—J)W = 2 | (DYDZ) dt, 


D(D?—I—J)W = 2421 | (YDZ) dt +27 | (ZDY) dt, 


x 


D*D*—I—J)W = 21’ [ (YDZ) dt +2J’ [ (ZDY) dt + 


+4/JW+([+J)(D?—I—J)wW, 
i.e., writing for brevity 
© = D?—I]—J, 
(O2—41J)W = 2I' [ (YDZ) dt +2J’ [ (ZDY) dt. 
On eliminating the integrals from (47)—(49) we have 
(I—J)(0?—41.J)W—(I’—J')DOW — 
—2(IJ'’—I'J)DW = —2(I'—J’). (50) 
Now operate with (J’—J’)D—(1”—J"): the right-hand side of (50) 
is eliminated and we have, as the required equation for W, 
[(1’—J’)D—(I"”—J") || (I—J)(0?— 41 J) — (I — J) DO — 
—2(IJ’—I'J)D|W = 0. (51) 
We note that, if the second bracket stood alone, we should have* 
the equation satisfied by the products of solutions of the equations 
of the second order only. 
* (2), 146. 
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The equation (51) may be more simply expressed as 
(I—J){{(I’—J") D—(I"— J") [0?— 41 J] W— 
—(I?—J") PW —2(I' J” —J' I") DW + (I—J)(I’— JPW} = 0, 
(52) 
a form admitting considerable further simplification in an important 
class of cases for which I’ J” = J’I". 


10. Whittaker’s functions 
As a final application of the methods of the earlier paragraphs we 
determine the equation satisfied by 


z= 9,(4;¢; px) Xy,(a';c; gx), 
where 
(4; ¢; px) = (px) exp(— i px) F(a; ¢; 2). 


Since y,(a;c; px) satisfies the differential equation 
d¢y = (p?+ka+ Px*)y, 
where k = }p(2a—c), P = 4p, p = }(c—1), we may, remembering 
that dy = dy/dw, where x= e”, derive the differential equation 
satisfied by z from the general result of §9. It is, however, somewhat 
less laborious to proceed as follows. 
If k’ = 4q(2a’—c), Q = 3g, then z satisfies the partial differential 
equations 
jz = (p?+ka+ P22x*)z, 
whence, by addition, 
(5°— 2p?)z—- (k-+-k’ )az—( 
Thus 
5(5°— 2p*)z—(k+-k’ )a(S+- 1)z—(P?+ Q?)x?(5+-2)z 
= 2466°2+ 20472, 


5(52—4p?)z—(k+k’ )a(5+ 1)z—(P?+ Q?)x?(5+ 2)z 
= Wz = 2(ka+ P22?)bz2+ 2(k’n+ Q*x*) Oz. 


2u(3—1){k8p-+k’S0}z-+ 22(8-+ 1){ P°84-+ Q780}z 
2ar(5—1){eb2-+-h'62-+ (kk )O$}2-+ 
+ 2x°(8+ 1){ P262-+ Q262-+-(P2+@Q2)0d}z. (56) 
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If now on the right we replace 6°z, ¢*z, 20¢z by their expressions 
taken from (54) and (55), then (56) is expressed in terms of the 
operator 6 only, and on rearrangement we obtain as the required 
equation 


8(5— 1)(8— 2)(8°— 4p?)z—(k-+-k’ )a 8(8— 1)(28+- 1)z— 
— 2(P?+ Q?)a? 8(5+ 1)?2+ (k—k’)?2?8?z++ 
+ (k—k’)( P?— Q?): (3(28-+-3)z+ 


If 
Y(a, Cc, px) = (px)-*¢-Yexp(— 4h px) ,F(a+1—c; 2—c; px), 


the general solution of (57) is 


x& 


’ |y,(a,c, par), Yo(a,c, px) Y,(a',c, gx), Y(a’,¢, qx) dt (58) 
'y,(a,¢, pt), yo(a,é, pt) || y,(a’,c, gt), yo(a’,c,qt)| t ! 


We may note that, employing the terminology of §1, the equation 
for y is a three-term equation while that for the product z has five 
terms. 
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INTEGRAL EQUATIONS FOR HEUN FUNCTIONS 
By A. ERDELYI (Edinburgh) 
[Received 4 March 1942] 


1. Some years ago C. G. Lambe and D. R. Ward (3) obtained 
integral equations for Heun polynomials, i.e. for polynomial solutions 
of the Fuchsian differential equation of the second order with four 
singularities. The present note contains an extension of the results 
of their paper. The new results include (i) integral equations for 
typical polynomial or transcendental solutions of the general Heun 
equation, i.e. with arbitrary values of the accessory parameter h, 
(ii) more general nuclei than those considered by Lambe and Ward, 
and (iii) the most general nucleus for Heun polynomials. The work 
presented here is analogous to and arose from certain investigations 
on periodic Lamé functions, which will be published in the Pro- 
ceedings of the Edinburgh Mathematical Society. Also there is a close 
connexion between the integral equations presented here and certain 
expansions of solutions of the Heun equation in series of hyper- 
geometric functions.* 


2. Let M, be the differential operator 


22 


M, = w(e—1)(e—a) + 


+ fy(e—1(a—a)-+82(¢—a)-+ex(e—1)}* +aBe, (1) 


and M, the differential operator obtained by replacing x by y. We 
assume throughout 

1+a+fp—y—i—e = 0. (2) 
The Heun equation is 

(M,,—aBh)u = 0. (3) 

Its typical solution which is analytical in the neighbourhood of and 
equal to unity at the origin is usually denoted by F(h,x), or, more 
fully, by F(a,h;«,8,y,8;2). Apart from exceptional (logarithmic) 
cases, the general solution of the Heun equation can be expressed 
in terms of functions of the type F(h, x), and therefore only integral 
equations for these functions will be considered. 


* For the latter expansions compare A. Erdélyi, Duke Math. Journal, 
9 (1942), 48-58. 
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The integral equations discussed will be of the form 


F(h,x)=2 ( y’(1—y)®1(1— y/a)* 1K (x, y) F(h, y) dy. (4) 
c 

It is seen either from the paper of Lambe and Ward or from the 
more general investigations of Ince (1), (2) that the nucleus K must 

satisfy the partial differential equation 
(M,—M,)K = 0, (5) 
and the contour of integration C must be such that the ‘integrated 
parts’ vanish. We notice that does not appear in (5); thus K may 
be taken to be independent of / and to depend solely on aq, f, y, 4, 


X,Y, @. . 
3. In order to find typical solutions of the partial differential 
equation (M,—M,,)K = 0, let us introduce new variables 6, ¢, where 


xy\} : .{(~—a)(y—a)) 3 
0s 6 = | —}*, s 08d = 7s 
cos (2) sin 6 cos ¢ ‘| naar | 


sin @sind = er. 1)\3 (6) 
| l—a 

Here 6 and ¢ may be conceived as polar coordinates on the unit 

sphere. Then (5) becomes 

eK 


sin?6 
: 002 


-+-{(1—2y)tan 6+ 2(5-+-«— })cot 6} a — 4K | 2: 
Cc 


eK , 1 

4 “~* 4 {(1—28)cot d—(1—2e)tan¢}— = 0. (7) 
og* | og 

This equation can be integrated by a product of a function of @ and 
a function of ¢. Typical solutions are 

0 ] OO ( 0 1 ° @) 
Pi 0 4-—8—c a cos%?@)xP; 0 0 —3+68+0 cos*d}, 

beg t—-eto B cas 1—8 —}+.e«-oc 

(8) 

where P{ } is the usual Riemannian symbol, o is an arbitrary 
separation constant, and 


xy 


cos?@ = —%, cos’é = ae (9) 


a 1—a)(xy—a)’ 
Two simple types of special nuclei will be mentioned. If co = 4—8, 


there is a branch of the second P-function which is constant; asso- 
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ciating it with that branch of the first P-function which is regular 
at the origin, the special nucleus 


F(a, B; y; xy/a) (10) 


is obtained. A similar simplification is possible for o = e—}. The 
nucleus (10) for « = 0, —1, —2,... is the nucleus discussed by Lambe 
and Ward. Again, if o is one of the roots of the equation, 


(5—}-+0—a)(8—}-+0—B)(e—}—o—a)(e—}—o—f) = 0, 


one branch of the first P-function in (8) is an elementary function, 

and we have nuclei which, apart from an elementary factor, are 
onan ie heen ; if : 1 

hypergeometric functions of («—a)(y—a)/{(1—a)(ay—a)}. 


4. Let us return now to the integral equation. Suppose that 
K(x,y) is a function of the form (8), or else a sum or integral with 
respect to o of such functions. Then 

[ y7(1—y)®4(1—y/a)e 1K (a, y) F(h, y) dy (11) 

é 
s a solution of the Heun equation provided that the ‘integrated 
parts’ vanish. This latter condition is satisfied, e.g., if C is a contour 
closed on the Riemannian surface of the integrand. Let us assume 
furthermore that the singularity y = a/x of the integrand is outside 
C, and that the branch regular at the origin of the first P-function 
in (8) has been chosen. Then (11) will represent an analytic function 
of x regular at the origin, where it takes the value 


[ yt —y)P4(.— ya) 7K 0, y) Fh, y) dy. (12) 


C 


Therefore (11) is a multiple of F(h, x), and (4) is true, 1/A being given 
by (12). It is important to remark that K(x,y) and C can always 
be chosen so that (12) does not vanish. In fact, a branch of the 
second P-function in (8) can be chosen so that 


(l—y)>1K (0, y) F(h, y) 


is analytic at y = 0 but not analytic at y = 1; with such a nucleus 
the Pochhammer double-loop slung round y = 0 and y=1 is a 
suitable contour unless y or 6 or both are positive integers: in the 
latter cases a simple loop or the interval (0,1) will be a suitable 
choice for C. 
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Thus F(h, x), with arbitrary values of h, satisfies the integral equa- 
tion (4) provided that K is either a nucleus of the form 


(xy—a))-S-4F (}—3—o-ba, 4—8—o+ By; 4 x 
a 


0 1 00 
xr @ 0 —}+6+o 
I—e_ 1—8 —}+e-—co 
or else a combination (sum or integral with respect to o) of such 
nuclei, and provided that C is a contour closed on the Riemann 
surface of the integrand, such that y = a/x lies outside it and that 
(12) does not vanish. 


(e—a)(y—a) 
(1—a)(«y—a) |’ 


(13) 


5. At first it may seem strange that we have obtained integral 
equations for F(h,2) without subjecting this function to boundary 
conditions and thereby selecting certain characteristic values of the 
accessory parameter h. The explanation is that (4) is, in general, 
a singular integral equation. Adjoining appropriate boundary con- 
ditions, it is possible to obtain from (4) integral equations of the 
usual Fredholm type. 

A typical pair of boundary conditions is regularity at x = 0 and 
at x =1. Solutions of (3) satisfying these boundary conditions will 
be called Heun functions, this name being analogous to that of 
(periodic) Lamé or Mathieu functions. Now, F(h,x) is already 
regular at « = 0. The condition that it should also be regular at 
x = 1 yields a transcendental equation for h. In the present section 
it will be assumed that / is one of the denumerable sequence of roots 
of this transcendental equation. It will also be assumed that in (13) 
the branch belonging to the exponent 0 at the singularity 1 of the 
P-function has been chosen. Under these assumptions the integrand 
of (11) has only multiplicative branch-points at y = 0 and y = 1. 
We choose C as a double-loop slung round 0 and 1, and, under the 
additional hypothesis that R(y) > 0, R(5) > 0, the contour integral 
will be a constant multiple of the integral extended from 0 to 1. 

Thus it is seen that, if (i) h satisfies the above-mentioned tran- 
scendental equation, (ii) R(y) > 0, R(S) > 0, and (iii) K is a nucleus 
of the form 

(ay—a)!->-° F(4—8—o+a, }—85—o-+f; y; 2y/a) x 


_s. Ux—1)(y—1) 
xF(—H +0, —bteoid ye). (s 
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or else a combination of such nuclei, then F(h, x) satisfies a Fredholm 
integral equation of the form (4) with 0 and 1 as limits. 

It is to be noted that no restriction has been imposed on a or f. 
Consequently, the solutions (eigen-functions) of the Fredholm integral 
equation are, in general, transcendental Heun functions, and not Heun 
polynomials. 


6. In this last section Heun polynomials will be discussed. Such 
polynomials arise if either « or f is an integer < 0. To fix ideas we 
shall assume throughout the rest of this note «a = —n where n is 
a non-negative integer. In this as in the general case all results of 
the previous section hold true without restriction or modification; the 
additional features will be discussed in this section. 

In this as in the general case there is a transcendental equation 
for h as a condition for the existence of Heun functions. The special 
feature in the case a = —n is that n+1 of the Heun functions are 
polynomials of degree n in x.* These polynomials are called Heun 
polynomials. Besides these polynomials—a fact very often over- 
looked—there is an infinity of transcendental Heun functions even 
for « = —n. 

Heun polynomials are naturally fully included in the general state- 
ment of the preceding section, and hence they satisfy integral equa- 
tions with nuclei of the type (14). In general, i.e. with arbitrary 
values of o, these integral equations will be satisfied by transcendental 
Heun functions as well as by polynomial ones. The question arises, 
and will be fully answered, whether there are nuclei all eigen- 
functions of which are Heun polynomials. 

From the theorem on the bilinear development of the symmetric 
nucleus of an integral equation it follows that, if there are such 
nuclei, they must be polynomials of degree n in x; and conversely, 
if K is a polynomial of degree n in 2, all solutions of (4) share the 
same property. There are exactly n+-1 such polynomial nuclei of 
form (14), corresponding to the n+1 values 


on, = 4—5—m (m=0,1 


(15) 


of a; o) corresponds to the nucleus discussed by Lambe and Ward, 
the other » nuclei are new. At the same time (15) exhausts all 
possibilities, or, more precisely, any nucleus all characteristic func- 


* Compare, for instance, (3) § 2.1. 
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tions of which are Heun polynomials is a linear combination, with 
constant coefficients, of the nuclei characterized by (15). The proof 
of this is immediate. 

This last theorem has interesting consequences. The product 
F (h, x)F, (h, y) of two Heun polynomials can be expressed as a linear 
combination of the n+1 special nuclei (15). Putting y = 0, 1, or a 
in this expression, three different expansions of a Heun polynomial 
in terms of Jacobi polynomials ensue. 

A great variety of related theorems may be obtained from the 
results of the present note by appropriate transformations of the 


variables in the Heun equation. 
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ON FUNCTIONS REPRESENTED BY 
CERTAIN SERIES 


By T. VIJAYARAGHAVAN (Dacca) 
[Received 21 July 1942] 
THEeorREM. Ifa, > 0, a, = O(1), a, 4 0 for an infinity of values of n, 


and 





n = an 
a, 2 a,2 


f(z) “i ? g(z) = > 
iil a all 
then the unit circle is the natural boundary of f(z) and also of g(z). 
This result is a sharpening of a theorem proved by G. Pélya.* 
That neither of the conditions a, > 0, a, = O(1) is entirely super- 
fluous can be seen from the following known examples: 




















p(n)2” ae... d(n)z” oe oes 22? 
1-2 (1—z)?’ 14-2" (1—z)? (1—2)” 
p(n)” ° p(n)z” =" 
ig? —_— —lUlCtC~ 


It may be remarked here that f(z) and g(z) are connected by the 
equation g(z) = f(z)—2f(z*), so that g(z) is a rational function of z 
whenever f(z) is a rational function of z. The converse is not true, 
for, if a, = n when n = 2” (m = 0, 1, 2,...) and a, = 0 otherwise, 
then it is easy to verify that | f(re?*742")| + 0 as r >1—0, k and m 
being any pair of integers; hence the unit circle is the natural 
boundary of f(z), whereas the corresponding g(z) is equal to z/(1—z). 
In proving the theorem of this note I concern myself with g(z) 
only, as the proof in the case of f(z) is similar and slightly easier too. 
It is easy to see that the series for g(z) is convergent for |z| < 1, and 
that its sum represents an analytic function of z. It is shown below 
that there exist arbitrarily large integers m such that every primitive 
(2m)th root of unity is a singular point of g(z). Actually it is proved 

that, if € is such a primitive root, then 
lim |g(ér)| = ©. (1) 

r—1—0 

From the prime number theorem (or from less deep results con- 
cerning the distribution of primes) we can deduce easily that, if « 
be any positive number and m any suitably large integer, then the 
primitive (2m)th roots of unity are so distributed on the unit circle 


* Proc. London Math. Soc. (2) 21 (1922), 36-8. 
3695.13 I 
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that there is at least one such root in every arc of length «. From 
this observation and (1) it will follow that the singular points of g(z) 
are dense everywhere on the unit circle, and hence that the unit 
circle is the natural boundary of g(z). The rest of this note is con- 
cerned with proving (1) for suitable values of €. 

Let g,,(€) denote the sum of those terms a, 2"/(1+-2") for which 
n is an odd multiple of m, and h,,(z) the sum of the remaining terms 
of the series for g(z). If € is a primitive (2m)th root of unity, 
0 <r <1, and v is not an odd multiple of m, then, by considering 
the position of the point (€r)" in the complex plane it is easy to 
verify that - 
\1+(Er)"| >c¢ = c(m) = 1—cos—, 
| 1 P 1 

1+(&r)"| ~c 


If M is the upper bound of a, a, a3,..., and 0 < r < 1, then we have 


and that therefore 


M< M 
ae SOS ae = (2) 
To estimate g,,(€r) we observe that, if n is an odd multiple of m, 


say n = m(2/-+-1), then 








1 aalEr™ _ (Lara ay cyanea MOY 
1+ -(ér)” l—-r”= ~*~ m(21+ 1) 
If now we suppose that 
. An(21+1) 
> a = 0, (3) 
9) 1 
SS m(2/-+-1) 
then, since a,, > 0, it follows that 
lim (r—1)g,,(ér) = 0. (4) 
r—>1-—0 


If (3) is true for an infinity of values of m, then (4) is true for those 
values of m, and then it follows from (2) that 


lim |9m(€r)| = lim |Jm(Er)+-Am(Er)| = 90 
r—+1—0 r—>1—0 


for a set of values of € which is dense everywhere on the unit circle. 
Now in order to complete the proof of the theorem it remains only 
to dispose of the case in which the expression on the left-hand side 
of (3) is convergent for all large values of m, i.e. 


S _Anfats)_ /— H+1. H4+2. H43....). 5 
aa (m +1, H+2, H+...) (5) 


i=0 
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It is proved below that in this case 
(1—r)h,(Er) > 0 (6) 
as r->1—0, m being any positive integer, and é being, as before, 


any primitive (2m)th root of unity. Moreover, 


lim (r—1)g,(ér) >“, (7) 
r—1-—0 m 


and, since a,, > 0 for an infinity of values of m, the theorem follows 
from (6) and (7); however, we have yet to prove (6), and for this 
purpose the following lemmas are needed. 


Lemma 1. If 


re) a © 
_— m(2l+1) < 0, e) a a ym(2l+1) 
then (1—r),,(r) > 0 as r >1—0. 


7 _ Gnlat+) __ 

For, J Yn(r) dr , $m Ipm(2i-+1) < < 5, < ©, 
and, since #,,(r) is a monotonic function of r, the conclusion of the 
lemma follows from the convergence of the integral. 

We define the density of a sequence ,, m., ”3,... to be the limit 
of m/n,,, provided that the limit exists as m— oo, and denote it by 
d(Ny, Ng, Ng,---)- 


Lemma 2. If «is any positive integer, and p,, Po,..., Pj, are odd num- 
bers which are prime to each other in pairs, and 0 < ny < Ny <3 < 
is precisely the set of integers that are such that none of them is an odd 
multiple of any number of the form 2“p, (0 << u<a;l<gu<k) then 


tmred=f-B piggy) 


Let s, denote the sequence of integers 
M-1p (214-1) (v =1, 2,..., &; t= 0, 1, 2,...) 


arranged in order of magnitude and without repetition. The se- 
quences 8}, 8»,..., 8, are all distinct (i.e. no two sequences have a 
common member) and, if d,, denotes the density of s,, then it is easy 


to see that 
sm al-(-alt-at-a 
2" Pp Pp Pr. 
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If we remove from the sequence of natural numbers the sequences 
81, Sg,..., 8,, then the resulting sequence is ”,. %, N3,... and its density 
is equal to 1—d,—d,—...—d,, which is plainly equal to the expres- 
sion on the right-hand side of (8). 


Lemma 3. If « and H are given positive numbers, then there exist 
Primes Py, Po,---, Py, and an integer « such that H < py < py < ... << Pp, 
and d(n,, so, 3,...) < €, where the sequence n1, No, Ns,... 18 defined as in 
the previous lemma. 

Let p, be a prime greater than H and pg, ps,... be the succeeding 


(1—F)(1-2}..(1-2) +0 
Pi Pe Pr 


as k > 00, it is seen from Lemma 2 that d(n,, 79, N3,...) < € Whenever 


primes. Since 


k and « are both suitably large. 


Lemma 4. If a, < M, 8 = d(ny,, no, N3,...) < €, and & is a primitive 
(2m)th root of unity, then 


lim (l—r) Sa, r" < Me. 
— nt 
rol—s 
If s, denotes > a,,, then it is plain that lims,/v < Me, and hence 


Nisv 


that 


(l—r) > a, 7" = (1—r)? 2 * r’° < Me 


p 


whenever 7 is sufficiently near to but less than unity. 


Lemma 5. (1—r)h,,(r&) > 0 as r > 1—0. 

Let 7 be any positive number. Determine +p, /y...., p, and « as 
in Lemma 3, taking H to be the same as in (5) and « = yc/M, where 
c¢ = 1—cos(z/m) and M is the upper bound of a,, dy, a@g,.... Since 
a, > 0, we have that (l—r) } a,7r" < },+},, where 

ie 9) 
ss = (1—r) 3 Pou, (7), } a = (1—r) Zz a,,7". 
0su<a t=1 
1sv<k 
From Lemma 1 we have that >, — 0 as r > 1—0, and from Lemmas 
3 and 4 we see that >, < Me = nc, asr >1—0. Hence, as r>1—0, 
a, 
\(1—r)hiy(ér)| <> Sern s 
Cc 
Since y is arbitrary, the lemma is proved and this completes the 
proof of the theorem also, 
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1. Tuts note contains a short statement of some new results con- 
cerning closed Hermitian transformations of deficiency-index (m, m)t 
and their resolvents in Hilbert space §. The most important of them 
seem to be those which are formulated in Theorems 5-7. They are 
derived from some new analytical representations stated in Theorems 
1-4. A more detailed development of the proofs will appear in a later 
paper, in which the case m = 00, omitted from the present note, will 
also be discussed. Detailed proofs of Theorems 6 and 7 are, however, 
given here because those in the later paper are of a different charac- 
ter. The abbreviations ‘c.H.t.’ and ‘d.i.’ will be used respectively 
for ‘closed Hermitian transformations’ and ‘deficiency-index’. 


I 
2. THEoREM 1. Let H be a c.H.t. of di. (m,m), H* its adjoint 


0 0 
transformation, H a self-adjoint extension of H, and R, = Hz", where 


0 0 
|. = H—xl, the resolvent of H. Let 4,(x), $o(x),..., b»(x) be a set of 
m linearly independent characteristic solutions of H* belonging to the 
characteristic value x, so that H* $,(%) = 0 (u = 1, 2,..., m) for every 
non-real x. 
Then the resolvent R, and the set $,(x) determine, for every non-real 
x, a matrix C(x) = {c,,(x)} of order m, which has an inverse matrix 
O-\(x) = {e,,\(a)}, 80 that 
m m 
(x—y)R, $,(x) = $,(%)— 2, 2, Curl® Cr (Y) PY) (1) 
v=l«K= 
for every non-real x and y. 
To prove this we first observe that there is a matrix A(z,y), 
whose elements depend on both variables x and y, such that 


(x—y) By by (0) = b4(2)— ¥ Ayal, Why). (2) 


+ The notation is that of Stone’s book (2). See (1) 87, definition 15; (2) 81, 
definition 2.21, and 338, definition 9.1. 
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We now apply the transformation R, to (2) and use the familiar 
fundamental equation for resolvents, viz. 

(z—y)R, R, = R,—R,. 
We thus obtain the matrix-equation 

A(x, y)A(y,z) = A(z, z). (3) 
If we substitute z = 7 in (3) and write C(x) for A(x,7), (3) becomes 
A(x, y) = C(x)C—(y), which is the desired result. 


3. If we put ®,, (x) 


formula (1) becomes 
(x—y)R, ®, (: ) _ ®,,(x)—®, (y). 

From this we immediately obtain 

THEOREM 2. The function {®,(x),9}, defined by (4) for any element 
g of §, is regular and analytic in every simply-connected domain of x 
which does not contain any point of the real axis, and 

a® ,(x)/dx = R,®, (x) 

is an element of § for every non-real x. 


0 
4. If H(A) denotes the resolution of the identity belonging to H, 
and AE the projection H(b)— E(a), where b > a, we obtain 


0 0 
H, AB®,,(x)—H, AE®,(y) = 0 


0 
by applying the transformation H, AE to (5). 
This formula leads to 


THEOREM 3. If we write 
b 


0 
YA = H,AE®,(x) = | (A—x)d{EQ)®,(x)}, (6) 
then ys), is an element of $ which is independent of x. 

If D denotes the domain of H, M, the linear manifold into which 
the projection AE takes §, then the yp (u = 1, 2,..., m) are elements of 
M,, and D.M, ts the set of all elements f of M, which are’ orthogonal 
to each of the ys. 


The last assertion in the theorem follows from the relation 


0 0 
0 = {H, AEf,®,(x)} = {AEH,f,®,(x)} = {f, H, AE®,(x)} = (f, 42), 
which holds for every element f = AEHf of D.M,. 
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5. THEOREM 4. 
m es 


jim |2 2 My wal = HY |A—2|? a{ BQ) Xu, 0 42) Pn) = 00 


for every system of m real or complex values of u,,. 


For suppose that this limit is finite for some system of numbers 


m . 0 
u,. Then ®(z) = 2 ta Pula) belongs to the domain of H, and so 


H,,®(x) = 0 for every non-real x. This is a contradiction, since the 


0 
self-adjoint transformation H cannot have any non-real charac- 
teristic value. 


6. We now divide the whole real axis into an infinity of intervals 
A, (—00 < p < +00) of finite length. Since®,(z)= S A, B®,(z), 
p=-@ 


we obtain from (6) an analytic representation of ®,,(x), namely 


< d{ E(A)ypAe 
®,, (x) = ee A, ER, Wr = S i d{ EA) Owe 4 
p=-@ 
This formula may also be written as 


7. Before passing to our three main theorems we need a new 
definition: a c.H.t. H of d.i. (m,m) will be called a prime transforma- 
tion if § does not contain any closed linear manifold which reduces 
H and with respect to which H is self-adjoint. 

In what follows we shall be concerned with prime transformations 
only. We can, however, easily reduce every c.H.t. H of d.i. (m,m) 
which is not prime to a prime transformation. Let I be the greatest 
closed linear manifold contained in § which reduces H and in which 
H is self-adjoint. Then H will be a prime transformation with respect 
to the Hilbert space §’ = § © M. 

Further, it is easy to show that, if H is a c.H. prime transformation 
of d.i. (m,m), and E(A) a resolution of the identity belonging to any 
self-adjoint extension of H, then the equation AZ®,(x) = 0, holding 
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for every » = 1, 2,..., m, for any non-real x, and for any interval A, 
implies AZg = 0 for every element g of §. 


8. THEoREM 5. In order that the spectrum of any self-adjoint exten- 
sion of a c.H. prime transformation H of d.i. (m,m) should contain 
a finite number of points only in every interval a+-8 < A < b—S, it is 
necessary and sufficient that the m elements $,(A+-in) should be con- 
tinuous in the open rectangle a < A <b, |n| <e. 

The fact that the condition is necessary follows from a theorem 
of Calkin} or from our formula (7). 

In order to prove that the condition is also sufficient, we first show 
that the continuity of the m elements ¢,(A+-7y) implies that of the 
m* functions c,,,(A+-77) in the same open rectangle. We then deduce, 
by means of Cauchy’s integral-representation, that, in the sense of 
Theorem 2, the m elements ®,,(x) are also analytic in the real points 
«x = X of the open interval a < A < b, except for a set of points A,, 
in which the determinant |/C(A,))|| vanishes, and which is at most 
denumerable. These points A, are poles of the ®,(), and are all the 


0 
characteristic values of H in the interval a < A < 6. We further 
notice that b 
° d{E(A)pa, i) 


_— pasag MDs 


{AE®, (x), y4} = | 


a 


is a regular analytic function for a << « < b, x #A,. This implies, 
however, that {H(A)pi, Wi} = |H(A)%7|\? remains constant in every 
interval A’ which does not contain any one of the roots A = 4, of 
||C(A)||. Hence A’E®, (x) = A’ Ey) = 0 for every x, and so we obtain 
finally A’E = 0, since H is supposed to be a prime transformation. 
Thus we have proved that the open interval a < A < b does not 


0 
contain any part of the continuous spectrum of H. 
0 

We prove finally that the number of characteristic values A, of H 
in every interval a+5 < A < b—&S is finite by considering the two 
mth-order determinants 

O(x, 2) = |{, (2), ®,(@)})|, $e) = ||fb (0), 4,8) 
We readily obtain the equation 
1 \|C(a)|| .||C(@) 


O(x,z) d(v,t) 
+ (3) 506, Theorem 3. 
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and then show that, within our open rectangle, 1/¢(7,Z) is a con- 
tinuous function of x = A+in and 1/®(x,Z) is an analytic function 
of x. Now it follows from (8) that, by a familiar elementary theorem 
of the theory of functions, ||C(zx)|| has only a finite number of roots 
in every interval a+65 < A < b—8: and this is the desired result. 

Theorem 5, of course, also yields the necessary and sufficient con- 
ditions that the spectrum of every self-adjoint extension of a c.H. 
prime transformation of d.i. (m,m) consists of its characteristic values 
only, which, moreover, have no finite limit point. 


Ill 


9. The existence of c.H. prime transformations of d.i. (m,m) whose 
self-adjoint extensions have spectra consisting of every point of the 
real axis becomes evident from the following two theorems. 


0 
THEOREM 6. Let H be a self-adjoint H.t. with simple spectrum, E(A) 


its resolution of the identity, and R,, its resolvent. Also let g be an 


0 
element of the domain of H such that the linear manifold M(J) spanned 
by E(Ag, if X runs from —oo to +00, coincides} with the whole Hilbert 
space $. Further, let p(A) = {B(A)g, g}. Then, in order to construct all 


0 
the c.H. prime transformations H of d.i. (1,1) to which H belongs as 
a self-adjoint extension,{ we must consider all those functions Q(A) 
(—o <A< +00) which fulfil the following three conditions: 


(i) the set of all the points A for which Q(A) = 0 is of p-measure zero; 


+0 


 f (QA)? eee ' 
(ii) | ¥o1 dp(A) converges; (iii) J |Q(A) |? dp(A) diverges. 


— 2 


Further, of we put 
7 2a 
oa) = | 2 aeony, (9) 


-— @ 


+ See (2) 243, Theorem 7.2, and 257, Theorem 7.9. 

t Von Neumann ((1) 108, Theorem 49) stated that every not-bounded H.t. 
T can be interpreted as the extension of another H.t. S. But he gave no 
method for the construction of all the transformations S having 7 as their 
extension, nor did he investigate the value of the deficiency-index of the 
transformation S constructed by him, which possibly might be infinite. 
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where x is any non-real number, we obtain the domain D which deter- 
mines the desired transformation H of d.i. (1,1) by applying the trans- 
formation R, to the linear manifold of all the elements of § which are 
orthogonal to ®(—7). We have finally 

H* Q(x) = 0 (10) 
for every non-real x. 

To show that these three conditions are necessary, we observe 
that formula (9) and condition (ii) follow, after a slight alteration, 
from the representation (7) for ®(x) by making use of the operational 
calculus; condition (iii) follows from Theorem 4, and condition (i) 
from the last paragraph of § 7. 


10. In order to prove that these three conditions are also sufficient, 
we have only to show that the domain D constructed in the way 
indicated in Theorem 6 is everywhere dense when these conditions 
are satisfied. 

If f is any element of § the operational calculus leads to the 
following representation of f by a function F(A): 


f= F(A) d{E(A)g}, (11) 


= [ IFO) dpQ). (12) 


On the other hand, any function F(A) whose integral (12) exists 
0 
furnishes an element of § by (11). Now we obtain the domain D of 


0 
H as the set of all elements 


f= Rf= | 5) aways (13) 


and the desired domain D as the set of all those elements 


T GA) 
‘ | i 


—@ 


d{ E(A)g} 


which satisfy the condition 


+o 


{h,®(—i)} = | ays 
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It is sufficient to show that the domain D is everywhere dense in 


0 0 
D, since D is itself everywhere dense in §. If we write 


O-) gy _; 
f + ODP (—4), 


where f is any element of § and {h,®(—i)} = 0, we obtain for every 


0 
element f’ of D the representation 


f= Rae HC ¥ R,®(—i) 


_ [ ® an 4.01— } 7 20) ene 
= f Maman LOC? | eS azo. cs) 


ae — 2 








Hence the only thing still requiring proof is that the element 


+0 


= | ee SO) d{BO9} (17) 


—-@ 





can be approximated to by an element (14) of D, since the other 
term of the expression (16) is already an element of D itself 

We first construct a complete ortho-normal set of functions Q,(A), 
Q,(A),... such that 


@ 


i: 20 ga) —s 





ry be (18) 


and O,(A) = QA{|O(i) |}. 


If, further, we put 








— ¥ 2,2,,(d) + 
G_.(a) = —2=3— ! K, a | ed 


G.,,(A) satisfies (15) because of (18). 
We now determine 25, %3,..., 2, in such a way that the expression 


ne} 


J Goi ad 


| R,®(-i) _,, 
{jo(—a PP ia] = 
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d{ E(A)g} 


a (19) 


is a minimum, which we denote by M,,. Since the M, evidently form 
a non-increasing sequence of positive numbers, lim M, exists. Thus 


n> 2 
the problem of approximating to the expression (17) is reduced to 


the question whether lim M,, = 


nD 


11. This, however, is answered by 


Lemma 1. Let (Ax,x) = y Ang Ep» X_ be a positive definite bounded 
pq=1 
n 
Hermitian form, (A™x,x) = Y ay,%,x, the nth reduced form, and 
pq=l1 
S(z) = {8,,(z)} = (A—zI)- the ete matrix. If M,, denotes the 


minimum of (Aa, 2) in the (n—1)-dimensional flat space x, = 1, then 


lim M, = 0 if and only if lim s,,(z) = 0, z taking negative values only. 
n yy “ie g 


n> x 


Proof. If S™(z) cae (AM—zI™) _— = {8 (z )} ‘ 
an elementary calculation leads to 1/M, = s(P(0). It is readily 
shown that for negative values of z, 
lim a)(z ) = 81,(2), an) (z j< aint P)(z ), 


n> 
sir(2) < sP2'), (2) < 8y(2’) 
eis : : ; 1 
and this implies lim M,, = lim - ; 
n> 20 834\2 


which is the desired result. 


12. In order to complete the proof of Theorem 6 we first notice 
that the expression (19) is the nth reduced form of a positive 
Hermitian form (Az,x) with x, = 1, 
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If we apply to the integral (20) the substitution » = ci and write 


Hap) = 2, |=) Hap) = 2( — J). 
nin foal PsA} t= fond A 


we obtain 


1 
nq _ | LH iprp(H)big(#) dw (1) + Pop(H)foq(H) darg(1)}. (21) 


Since the Q,(A) are supposed to form a complete ortho-normal 
system with respect to dp(A)/(A?+-1), it is readily seen that (21) gives 
the spectral representation of the Hermitian form (Az, x). 

This implies further that the resolvent s,,,(z) of the Hermitian form 
are defined by (20) or (21) is given by the 


whose coefficients a,, 


equation 
1 


8pq(2) = [ gan al) daw, (u) + Yop (1 )pog() darg(u)}. 


0 ' 
If we again substitute » = (A?+-1)-1, we obtain 
ay — [Mp AQ) 
Snq(2) _ [ 1—2(2-+1) dp(A). 
Hence we see that the condition 
PF QA)? 
ineti— tie | eee #am 


= @ 


and condition (iii) of Theorem 6 coincide. 
Finally, the assertion (10) follows from the equation 


{H;h’,®(x)} = 0, 
which holds for every h’ of D and can readily be derived from (9). 
This completes the proof of Theorem 6. 


0 
13. THrorEM 7. Let H, R,, E(A), 9, and p(A) have the same meaning 
as in Theorem 6. Then, in order to construct all the c.H. prime trans- 


0 
formations H of d.i. (m,m) to which H belongs as a self-adjoint 
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extension, we must rat all the sets of m linearly independent func- 
tions (2,,(A) (u = 1, m) which fulfil the following three conditions: 


(i) the set of all points in which the system of m equations Q,,(A) = 0 
is satisfied 1s of p-measure zero; 


*|Q, (A) |? 
v1 


—@ 


(ii) the m sacle dp(A) converge; 


(iii) the integral 9 | Su, dp(A) diverges for every system of 


m real or complex aatan of 7] 
Further, if we put 


Te 


= d{E(N)9} (uw = 1, 2,..-, m), (23) 


where x is any non-real number, we obtain the domain D which deter- 
mines the desired transformation H of di. (m,m) by applying the 
transformation R, to the linear manifold of all the elements of $ which 
are orthogonal to each of the m elements D,(—i). We have finally 
AF®,,(x) = 0 for every u and every non-real x. 


14. The proof follows the same lines as that of Theorem 6. 

We first orthogonalize the m elements ®,(7) defined in (23) by 
E. Schmidt’s familiar procedure in such a way that the equations 
(18) are satisfied for 1< p, g< _m. With this the linear manifold 
spanned by the m elements ®,(7) remains unaltered. We further 
construct the other functions Q,(A) (p = m+1, m+2,...) so that 
equations (18) hold for every p and q, and so that we obtain a com- 
dp(A) 
+1" 

The same argument as in §10 leads to the expressions correspond- 
ing to (19), namely 


plete ortho-normal set of functions Q,,(A) with respect to 


«oO 


v 2,0)+ x, 12,0) 
= |e ___ p= 2. a ee 
2) —h,|? = |: (A2-+- 1)? dp(A) 
(wu = 1, 2,...,m; nm >m-+1) 


whose minima J/# are now to be considered. 
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15. The question whether lim M+ = 0 is answered by 


Lemma 2. Let (Ax,x), (A™x,x), and S(z) have the same meaning 
as in Lemma 1. If, for n >m+1, M* denotes the minimum of 
(A™x, x) in the (n—m)-dimensional flat space 

x, =1, %=—%=—... = 2, = 0, (24) 
then lim M* = 0 if 


n—-> 2 


lim >) 8,,(z)t,uU, = 00 (z< 0) (25) 


2-0 pv=1 
for every system of complete numbers u,, for which 


2 Me han, (26) 


Proof. Let M,,(u,, U,..-,>U%,) = M,(u) denote the minimum of 
(Az, x) in the (n—1)-dimensional flat space 


Set =], (27) 


where the u,, satisfy (26). Then Lemma | implies that lim M,(u) = 0 
n> 


if and only if (25) is satisfied. This can readily be shown by sub- 
jecting (A™z,2) to a unitary transformation which carries the flat 
space (27) into z, = 1. We now consider the maximum yp, of M,,(u) 
on the unit sphere (26) and prove by a familiar argument that lim p, 


no 


is also zero if (25) is satisfied. This implies that 

M,(u)<e« forn > N(e) (28) 
uniformly for any uw on the unit sphere (26). 

If M* has the meaning defined in Lemma 2, the ellipsoid 
(A™z,x) = M*® 

touches the flat space (24) in the point at which the minimum is 
attained. Hence the (n—1)-dimensional tangent plane to the ellip- 
soid at this point must contain the (n—m)-dimensional flat space 
(24). Hence its equation has the form 


a ty, = Uy, (29) 


- where the w, satisfy (26). This implies, however, that VM is also the 
minimum of (Az, x) in the (n— 1)-dimensional flat space (29). More- 
over, since |u,| < 1, because of (26), we have 


Mx < M,(u) 
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for the special system of values wu, defined by (29). From this and 
(28) we obtain lim M@* = 0, the desired result. 


n—> x 
16. In order to complete the proof of Theorem 7 we again have 
to consider the resolvent S(z) = s,,,(z) of the Hermitian form (Az, x) 
whose coefficients are defined by (20). Then the s,,,(z) are determined 
by (22), and we obtain 
' m 9 
x sy ~ 
2X My Q2y() 
oh a-Aeh 
1—2(A2+1) 
Hence we see that conditions (25) and (iii) of Theorem 7 coincide. 
Finally, the assertion H¥ ®,,(x) = 0 (u = 1, 2,...,m) may be proved 
by the argument used in the last paragraph of §12. This completes 
the proof of Theorem 7. 
17. It may be remarked that a theorem similar to Theorems 6 


0 
and 7 also holds when the spectrum of H is not simple but of multi- 
plicity k, where k is a finite integer. We notice, however, that in 
this case we always have m > k for every c.H. prime transformation 


0 
H of d.i. (m,m) having H as one of its self-adjoint extensions. 


18. The application of the above methods to some special cases 
will be dealt with in a later note. This will give explicit expressions 
for the resolvents of all the self-adjoint extensions of a matrix J™ 
of d.i. (m,m) which is composed in a suitable way of m Jacobit 
matrices J,,, each of d.i. (1,1). In this case the hypothesis of Theorem 
5 will be satisfied. 

+ See (2) 530-614, especially 545, Theorem 10.27. 
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